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1 .  INTRODUCTION 

Reentry  vehicle  erosion  is  an  adverse  environmental  impact  on  weapons  delivery 
systems  with  potential  significance.  This  phenomenon  which  directly  affects  the  accuracy 
of  reentry  trajectories  is  attributable  to  mechanical  ablation  of  vehicle  aerodynamic  surfaces 
during  reentry  due  to  interaction  with  atmospheric  ice  clouds  and  solid  or  liquid 
hydrometeors  (i.e.  precipitation).  Assessment  of  existing  or  prediction  of  future  erosion 
severity  potential  requires  an  understanding  of  the  temporal  and  spatial  distribution  (i.e. 
climatology)  of  contributing  meteorological  conditions.  'ITiese  include  both  high  altitude 
cirnis  cloud  and  convective  activity  at  lower  altitudes.  'Ihus,  appropriate  analyses  and 
observation  techniques  are  necessary  to  characterize  cloud  and  hydrometeors.  This  report 
focuses  on  the  latter  issue.  In  particular,  due  to  the  site  specific  nature  of  the  erosion 
problem,  analyses  are  required  at  specific  land  based  locations. 

Significant  effort  has  been  devoted  to  the  development  of  analysis  techniques  based 
on  time/altitude  cross  section  analyses  to  identify  regions  of  likely  convective  .activity 
(Feteris  et  al.,  1976;  Hardy,  1979).  'Ihese  approaches  supported  by  radiosonde  and 
surface  station  data  were  labor  intensive  and  subject  to  the  sparcity  of  available  surface  and 
upper  air  data.  For  the  latter  reason,  for  example,  questions  of  spatial  representativeness 
are  always  present  (Bunting  and  Touart,  1980).  Due  to  the  level  of  effort  involved  in 
implementing  these  approaches,  they  are  unsuitable  for  application  to  large  scale 
climatological  studies. 

Both  cloud  and  hydrometeors  can  be  observed  globally  using  satellite  based  sensors 
(Isaacs  et  al.,  1986).  Previous  application  of  visible  and  infrared  satellite  imagery  data  to 
the  cloud  problem  is  discussed  by  Conover  and  Bunting  (1977).  The  advent  of  remotely 
sensed  microwave  imagery  such  as  that  from  die  Defense  Meteorological  Satellite  Program 
(DMSP)  Special  Sensor  Microwave/Imager  (SSM/I)  provides  die  capability  to  directly 
sense  hydroineteor  liquid  water  content  (Savage  et  al.,  1987).  Due  to  the  radiative  transfer 
physics  relevant  to  the  remote  sensing  of  precipitation,  the  phenomenology  of  rain  rate 
retrieval  is  different  over  land  and  oceans.  This  is  due  to  the  inherent  differences  in  land 
type  and  ocean  surface  emissivities  (Isaacs  et  al.,  1989).  Retrieval  techniques  have  been 
identified  to  provide  quantitative  measurements  of  rainfall  rates  over  the  ocean  based  on 
microwa  ve  imager  data  (Wilheit  cl  al.,  1977)  and  oceanic  climatologies  have  been 
developed  (NASA,  1976). 

Due  to  the  spatial  variation  of  emissivities  over  land  and  temporal  features  such  as 
snow  (which  have  signatures  similar  to  precipitation),  there  are  caveats  related  to  the 
retrieval  of  rain  rate  over  land.  However,  in  optically  thick  situations  when  surface 
contributions  to  measured  brightness  temperatures  are  small,  precipitation  can  be  inferred. 
Such  optically  thick  cases  are  those  which  characterize  intense  convection  and  these  can  be 
mapped  over  land  (Spencer  and  Santek,  1 985).  The  recently  completed  SSM/I 
calibration/validation  study  provides  further  confidence  in  the  validity  of  seeking 
quantitative  inferences  of  rain  rate  over  both  ocean  and  land.  The  specific  analysis  tool 
available  in  this  regard  is  the  rate  rate  retrieval  algorithm  developed  by  the  University  of 
Wisconsin  research  group  (Olson  et  al.,  1988). 

While  rainfall  rate  brightness  temperature  relationships  have  been  validated  based 
on  the  MJCS  154-86  requirements  for  precipitation  measurements,  the  physics  of  the 
reentry  erosion  phenomenon  indicates  that  a  knowledge  of  die  liquid  water  content  profile 
is  necessary.  This  can  be  achieved  by  adoption  of  a  suitable  model  of  the  vertical 
distribution  of  liquid  water  content  resulting  in  the  desired  surface  rainfall  rate.  'Hie 
selection  of  the  meteorologically  appropriate  liquid  water  content  vertical  distribution  model 
is  important  due  to  the  relationship  between  liquid  wafer  content  and  radiometric  signature 
through  the  dropsize  distribution  (Falcone  et  al.,  1979).  For  example,  many  investigators 
commonly  use  the  Marshal!  Palmer  (1948)  relationship  which  was  derived  for  stratifonn 
rain.  Erosion  severity  is  a  function  of  precipitation  severity  and,  therefore,  convective 
situations  are  of  greater  potential  interest.  For  these  cases,  other  relationships  are  likely 
valid  (Ulbrich,  1983;  Willis  and  Tattleinan,  1989). 
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This  research  study  applies  SSM/I  based  precipitation  r;mote  sensing  technology  to 
the  characterization  of  reentry  vehicle  erosion  environments  over  selected  sites.  Three 
complementary  areas  of  investigation  were  identified  in  our  study  proposal:  (a)  data 
acquisition  and  application  of  tire  SSM/I  rain  rate  retrieval  algorithm  to  the  development  of  a 
four  month  summer  climatology  of  rain  rate  over  selected  sites  with  an  emphasis  on  tire 
identification  and  analysis  of  heavy  rain  rate  situations,  (b)  identification  and  assessment  of 
appropriate  convective  cloud  models  to  establish  a  parameterization  of  the  relationship 
between  surface  rainfall  rate  and  the  vertical  profile  of  hydrometeor  liquid  water  content, 
and  (c)  testing  of  the  parameterization  of  hydrometeor  liquid  water  content  on  the 
determination  of  liquid  water  content  profiles  for  selected  ca  es  with  emphasis  on  intense 
convection.  These  issues  are  explored  m  the  following  sections. 

2.  BACKGROUND 

2. 1  Rain  Rate  Retrieval  Concepts 

At  microwave  wavelengths,  precipitation  sized  droplets  provide  a  source  of 
atmospheric  attenuation  analogous  to  the  effect  of  cloud  droplets  in  the  infrared  spectrum 
(Savage,  1978;  Falcone  et  al.,  1979).  This  attenuation  mechanism  suggests  a  direct  causal 
relationship  between  rainfall  and  microwave  atmospheric  emission  which  has  been 
exploited  to  infer  rainfall  rate  (Weinman  and  Wilheit,  1981).  Furthermore,  the  presence  of 
precipitation  within  the  field  of  view  of  mk'owave  sensors  is  itself  of  considerable  interest 
since  the  quality  of  resultant  retrievals  of  other  quantities  is  most  certainly  degraded  (Liou  et 
al.,  1981).  For  example,  microwave  sounding  brightness  temperatures  such  as  those  from 
the  SSM/T  must  be  corrected  for  rain  in  much  the  same  way  as  infrared  radiances  are  for 
cloud.  To  provide  a  theoretical  microwave  brightness  temperature/rainfall  rate  relationship 
applicable  to  rain  rate  retrieval,  models  of  both  the  rain  layer  (including  a  rain  layer  height 
and  thickness)  and  of  rain  microphysics  are  required.  These  models  determine  the  vertical 
distribution  of  hydrometeor  liquid  water  content.  Both  of  these  factors,  of  course,  vary 
synopticaJly  introducing  some  uncertainty  into  the  retrieval  process. 

Collectively,  these  factors  determine  the  atmosphere’s  scattering  and  absorption 
contribution  to  total  satellite  incident  brightness  temperature  which  is  evaluated  as  a  multiple 
scattering  calculation  assuming  that  the  rain  layer  fills  the  field  of  view.  Microwave 
radiative  transfer  theory  accounting  for  multiple  scattering  has  been  extensively  applied  to 
the  study  of  atmospheric  precipitation  and  clouds  (Weinman  and  Guetter,  1977;  Wilheit  et 
al.,  1977;  Tsang  and  Kong,  1977;  Jin  and  Isaacs,  1985).  The  surface  background  against 
which  the  rain  is  modeled  consists  of  surface  emission  and  surface  reflected  atmospheric 
contributions,  which  must  be  included. 

Over  the  ocean,  surface  emissivity  is  low  and  relatively  uniform,  providing  good 
contrast  for  the  quantification  of  precipitation.  This  approach  was  applied  over  the  ocean 
by  Wilheit  et  al.  (1977)  to  data  from  the  NIMBUS  5  Electronically  Scanning  Microwave 
Spectrometer  (ESMR)  operating  at  19.35  GHz.  Accuracies  of  2  mm/h  and  a  dynamic 
range  up  to  about  20  mni/h  were  achieved  using  ground  based  radar  for  validation.  Results 
based  on  retrievals  obtained  using  the  Seasat  and  Nimbus  7  SMMR  instruments,  for 
example,  indicate  that  the  retrieved  rainfall  rates  generally  underestimate  those  measured  at 
the  surface  with  rain  gages  (cf.  Gloerson  et  al.,  1984).  Lipes  (1982)  found  that  the  Seasat 
SMMR  failed  to  detect  showery  precipitation  associated  with  convective  cloud  in 
midlatitudes  and  underestimated  rain  rates  for  heavy  precipitation.  This  behavior  was 
attributed  to  both  loss  of  incremental  sensitivity  to  precipitation  at  higher  rainfall  rates  and 
insufficient  sensor  resolution  (about  30  km  at  37  GHz)  in  convective  situations.  Sensor 
resolution  plays  a  role  since  intense  precipitating  cells  with  horizontal  scales  of  a  few 
kilometers  will  not  fill  the  microwave  radiometer's  field  of  view.  Similar  results  were 
obtained  in  midlatitudes  by  Alishouse  (1983). 

Theoretical  calculations  suggest  that  overland,  higher  frequency,  dual  polarized 
measurements  can  distinguish  atmospheric  scattering  due  to  rain  from  surface  contributions 
(Savage  and  Weinman,  1975;  Weinman  and  Guetter,  1977;  Huang  and  Liou,  1983). 
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I  lowever,  in  practice  the  dynamic  range  of  measurable  rainfall  rates  is  so  much  reduced  at 
37  GHz  that  these  measurements  are  virtually  useless  and  considerable  ambiguity  of 
surface  and  atmospheric  contributions  still  exists.  In  addition  to  tire  problems  of  specifying 
appropriate  rain  models  and  obtaining  measurements  over  land,  another  problem  arises 
because  the  typical  horizontal  scale  of  precipitating  cloud  elements  is  generally  much  less 
than  die  field  of  view  of  the  ESMR  instrument  (50  km).  Integration  over  a  large  field  of 
view  creates  a  negative  bias  which  underestimates  instantaneous  rainfall  rates.  In  spite  of 
these  difficulties,  areas  of  intense  precipitation  can  be  identified  (Spencer  and  Santek, 

1985).  Based  upon  the  above  discussion,  there  will  be  some  difficulties  in  providing 
accurate  rain  rates  in  intense  convective  situations  from  microwave  data  alone.  The  rainfall 
rates  in  the  higher  ranges  of  the  desired  domain  will  almost  always  be  due  to  energetic 
convective  cells  with  spatial  scales  on  the  order  of  a  few  kilometers  which  cannot  be  fully 
resolved  with  the  25  km  FOVs  of  the  SSM/I  instrument.  In  such  cases  it  has  been  noted 
that  high  rain  rates  are  always  associated  with  minima  in  the  equivalent  blackbody 
temperature  (EBBT)  field  observed  from  infrared  imagers  (Negri  and  Adler,  1987a,b). 

'Hie  converse  of  this  is  not  true,  however,  i.e,  low  EBBTs  can  be  due  to  other  than 
precipitation,  therefore,  it  is  advantageous  to  exploit  both  microwave  (for  lower  rain  rates 
and  to  identify  the  presence  of  intense  convective  precipitation)  and  infrared  (at  higher  rain 
rates)  data.  At  the  higher  rain  rates  (and  to  help  in  the  determination  of  beam  filling  at 
lower  rates),  OLS  infrared  imager  data  can  be  used.  The  specific  approach  adopted  could 
be  based  on  that  described  by  Adler  and  Negri,  1988.  They  used  GOES  imager  data  to 
delineate  convective  rain  areas  by  searching  for  minima  in  the  EBBT  field  and  then 
assigned  rain  rates  based  on  the  results  from  a  one  dimensional  cloud  model  which 
provided  the  relationship  between  convective  development  (cloud  top  height)  and  the 
resulting  rain  rate.  Stratiform  rain  rates  were  delineated  based  on  EBBT  threshold  criteria. 
Due  to  level  of  effort  constraints,  however,  this  study  will  focus  on  die  use  of  SSM/I  data 
alone  to  derive  convective  rainfall  climatology. 

2.2  SSM/I  Data  Characteristics 

The  Defense  Meteorological  Satellite  Program  (DMSP)  Special  Sensor 
Microwave/Imager  (SSM/I),  launched  on  19  June  1987,  is  a  passive  microwave  radiometer 
which  provides  brightness  temperature  data  of  particular  relevance  to  the  monitoring  of 
global  precipitation  properties.  The  SSM/I  sensor  antenna  observes  seven  microwave 
channels  (three  dual  polarized  frequencies:  19.35,  37.0,  and  85.5  GHz  and  a  vertically 
polarized  channel  at  22.235  GHz)  and  scans  conically  with  an  angle  of  incidence  on  die 
surface  of  the  earth  of  53.1  degrees  (Savage  et  al.,  1987).  Using  a  single  antenna  to  span 
this  frequency  range  results  in  a  sensor  field-of-view  (FOV)  which  varies  with  frequency: 
about  50  km  at  the  two  lowest  frequencies,  25  km  at  37  GHz,  and  13  km  at  85.5  GHz. 

The  exact  spatial  resolution  of  each  FOV  depends  on  die  definition  used  and  the 
antenna  response  pattern.  Based  on  the  polar  orbiting  DMSP  satellite,  SSM/I  data 
potentially  provides  global  coverage  daily,  with  twice  daily  coverage  possible  at  northern 
latitudes  due  to  orbital  overlap.  SSM/I  data  are  archived  through  the  DoD-NOAA  Shared 
Metsat  Processing  agreement  by  NOAA/  NESDIS  through  the  Cryospheric  Data 
Management  System  (CDMS)  at  the  National  Snow  and  Ice  Data  Center  (see  Weaver  et  al., 
1987).  These  data  consist  of  bodi  satellite  data  records  (SDRs,  i.e.  the  calibrated 
brightness  temperatures)  and  environmental  d  a  records  (EDRs,  i.e.  the  retrieved 
parameters).  Table  1  illustrates  the  sensitivity  SSM/I  microwave  imager  channels  to  a 
variety  of  desired  geophysical  parameters  (the  DMSP  assigned  priorities  for  these 
parameters  are  in  brackets). 

The  retrieval  of  precipitation  from  SSM/I  data  foiiows  die  statistical  method  outlined 
in  Lo  (1983)  and  used  in  simulation  for  SSM/I  data  sets  in  Jin  and  Isaacs  (1987).  While 
simulation  retrieval  results  suggest  that  the  required  accuracies  of  a  few  mm/h  are  possible 
at  the  lower  rain  rates  a  number  of  important  factors  are  often  neglected.  These  include  the 
degree  of  beamfilling,  the  vertical  extent  of  the  precipitation,  the  rain  drop  size  distribution, 
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Table  1 .  Microwave  Iinager  Channel  Applications 


Sensor  Parameters 

Frequency  (GHz)  19.35 

Spatial  resolution  (km)  50 

Sensitivity  (K)  0.6 

Geophysical  Parameters, 

[PriorityJ 

Precipitation  (land)  [5J  o 

Precipitation  (ocean)  15]  • 

Snow  cover  [10J  • 

Sea  ice  (extent,  type)  [  1.8]  • 

Sea  surface  temperature  [8]  o 

W  ind  speed  (ocean)  [4 1  o 

Atmospheric  water  (total)  [3]  • 

Soil  Moisture  [9]  o 

Vegetation  (e.g.  Aibedo  [1 J  x 

Cloud  Liquid  Water  |7]  x 


22.235 

37 

85.5 

50 

25 

12.5 

0.8 

0.8 

1.0 

• 

• 

0 

X 

X 

0 

X 

• 

• 

0 

X 

X 

X 

X 

• 

0 

0 

X 

X 

X 

X 

• 

Channels  which  are  important  for  determining  each  parameter  are  indicated  using  the  following  code: 
♦  =  Critical;  O  =  Important;  x  =  Helpful. 


the  emperature  profile  through  the  precipitating  layer,  and  the  presence  of  glaciated 
precipitation  which  can  significantly  alter  the  brightness  temperature  signature.  These 
factors  are  often  kept  constant  in  simulations  whereas  they  vary  in  the  natural  atmosphere. 
Uncertainties  associated  wit!,  these  factors  certainly  impact  the  accuracy  of  rainfall  rate 
determination,  and  it  is  perhaps  more  reasonable  to  expect  that  a  few  broad  rainfall  rate 
categories  are  retrievable  from  the  SSM/I.  Retrievals  of  precipitation  (both  liquid  and 
glaciated)  and  surface  emissivity  from  simulated  SSM/I  data  were  discussed  in  a  journal 
article  by  Jin  and  Isaacs  (1987)  which  also  described  a  specially  developed  multiple 
scattering  mode!  which  was  designed  to  simulate  dual  polarized  brightness  temperatures  in 
the  presence  of  inhomogeneous,  nonisothermal  distributions  of  atmospheric  precipitation. 

The  capability  to  simulate  inhomogeneous  (i.e,  those  varying  with  height) 
distributions  of  precipitation  is  necessary  to  treat  the  realistic  variation  of  rain  rate  with 
altitude  within  developing  frontal  systems  and  the  phase  change  (from  water  to  ice) 
occurring  in  convective  situations.  In  that  paper  brightness  temperature  simulations 
specifically  applicable  to  the  SSM/I  were  shown.  Figure  1 ,  for  example,  illustrates  the 
dependence  of  simulated  dual  polarized  85.5  GHz  brightness  temperature  on  rain  rate  and 
the  presence  of  an  upper  layer  of  frozen  precipitation.  The  figure  inset  illustrates  a  model 
of  the  simulated  atmosphere  with  a  5  km  rain  layer  over  the  ocetui  surface  and  either  a  3  km 
ice  layer  or  another  3  km  rain  layer  above.  It  can  be  seen  from  these  results  that  due  to 
enhanced  multiple  scattering  at  this  frequency,  the  brightness  temperature  is  significantly 
lowered  by  the  ice  layer.  At  lower  frequencies  such  as  19.35  GHz,  the  ice  layer  has  little 
or  no  effect  and,  therefore,  using  the  multispectral  SSM/I  data,  the  phase  of  the  upper 
levels  of  precipitation  can  be  identified  in  the  retrieval/analysis  procedure  in  addition  to  the 
rain  rate.  This  provides  a  method  to  probe  the  vertical  structure  of  the  precipitation. 

A  significant  calibration/validation  effort  (Olson  et  al.,  1988)  has  focused  on  the 
SSM/I  precipitation  retrieval  algorithm.  The  resulting  investigation  provides  simple 
formulae  to  obtain  rainfall  rates  over  ocean  and  land  and  addresses  the  loss  of  the  85.5  V 
channel  data.  The  latter  issue  is  of  significance  for  the  inference  of  rain  rate  over  land.  We 
propose  to  employ  the  most  recent  land  regression  coefficients  obtained  from  the 
University  of  Wisconsin  group.  The  steps  in  the  retrieval  are  essentially:  ( 1 )  antenna 
patent  correct,  (2)  apply  existing  SSM/I  precipitation  screening  logic,  arid  (3)  apply 
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Figure  l .  SSM/I  dual  polarized  85.5  GHz  brightness  temperature  vs.  rain  rate  (min/lir) 
with  and  without  top  layer  of  glaciated  precipitation 

regression  formulae  for  rain  rate.  Step  (1)  is  part  of  the  coding  available  with  the 
compacted  SSM/I  data  from  Remote  Sensing  System,  Inc.  (F.  Wentz,  personal 
communication).  The  coding  for  steps  2  and  3  will  be  obtained  via  GL  from  the  University 
of  Wisconsin.  We  note  that  due  to  the  additional  computational  load  imposed  by  the 
antenna  pattern  correction  step,  it  might  be  prudent  to  develop  a  prescreening  algorithm  to 
identify  precipitating  situations  directly  from  the  antenna  temperatures  themselves.  If  this  is 
possible,  only  those  cases  identified  need  be  antenna  pattern  corrected. 

2.3  Convective  Cloud  Models 

As  discussed  in  the  previous  section,  the  SSM/I  retrieval  algorithm  provides  rain 
rate  as  an  output  parameter  although  the  radiative  transfer  physics  indicates  that  radiometer 
brightness  temperatures  are  fundamentally  sensitive  to  hydrometer  liquid  water  content 
(LWC).  Assessment  of  reentry  vehicle  erosion  indices  require  determination  of  a 
hydrometeor  liquid  water  content  vertical  profile.  The  required  interface  between  surface 
rainfall  rate  and  LWC  profile  is  a  model  of  the  vertical  distribution  of  hydrometeors 
appropriate  to  the  observed  synoptic  situation.  Previous  studies  such  as  Peirce  et  al. 

(1975)  have  focused  on  the  applicability  of  specific  hydrometer  LWC  models  for  this 
purpose.  Notably,  Falcone  et  al.  (1979)  have  specified  climatological  cloud  and 
hydrometeor  liquid  water  content  models  applicable  to  the  simulation  of  microwave  and 
millimeter  wave  data  sets  which  provide  strawman  cloud  model  candidates  as  functions  of 
precipitation  intensity  categories  (i.e.  light,  moderate,  heavy).  Since  the  application  driven 
climatological  emphasis  will  be  on  the  delineation  of  convective  precipitation,  our  focus 
will  be  on  convective  cloud  models. 

Table  2  provides  a  list  of  candidate  models  which  characterize  a  range  of  synoptic 
situations,  including:  frontal  rain,  thunderstorms,  tropical  squall  lines,  and  cumulus 
towers.  For  example,  a  simple  one-dimensional  convective  cloud  model  is  that  based  on 
the  work  of  Cotton  (1972a,b)  and  Simpson  and  Wiggert  (1969).  The  model  dynamics  are 
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Table  2.  Candidate  Convective  Cloud  Models 


Synoptic  Situation 
Thunderstorm 
Frontal  Rainband 
Cumulus  'lower 
Tropical  Squall  Line 
Tropical  Storm 
Convective  Cloud 


Reference 

Heltnsfield  and  Fulton,  1988 
Rutledge  and  Hobbs,  1984 
Simpson  and  Wiggert,  1969 
Tao  and  Simpson,  1989 
Wilheit  et  al.,  1982 
Cotton,  1972a,b 


based  on  the  conservation  of  vertical  momentum  including  the  buoyancy  effects  of 
perturbation  temperature  and  liquid  water  loading.  Physical  processes  included  in  this 
simple  model  are  the  latent  heating  by  condensation,  entrainment  of  environmental  air, 
conversion  of  cloud  liquid  water  to  rain  (and  its  fallout),  and  the  freezing  of  condensate. 
The  microphysical  parameterizations  are  kept  intentionally  simple,  both  in  the  interests  of 
run  time  and  because  of  the  likely  uncertainties  in  the  input  data.  In  Section  S,  we  discuss 
the  hydrometeor  cloud  model  in  greater  detail  and  identify  parameterization  for  stratiform 
and  convective  precipitation. 

3.  TECHNICAL  APPROACH 

3.1  Rain  Rate  Estimation  Procedure 

A  functional  flow  diagram  of  the  project  analysis  approach  is  shown  in  Figure  2. 
The  essential  steps  are:  (a)  reading  of  the  SSM/I  TDRs  from  the  acquired  SSM/I  data 
tapes,  (b)  calibration  and  antenna  pattern  corrections  applied  to  the  data  to  obtain  brightness 


Figure  2.  Functional  flow  diagram  for  SSMII  data  based 
erosion  parameter  estimation  procedure 


\ 
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temperatures,  (c)  geographical  binning  of  the  data  according  to  the  location  or  each  desired 
site  by  latitude  and  longitude,  (d)  application  of  the  rain  rate  retrieval  algorithm  to  the  data, 
and  (e)  evaluation  of  relevant  climatologies.  The  last  step  includes  calculation  of  time  series 
for  the  four  month  period  for  the  rain  rate  and  integrated  liquid  water  content  of 
precipitation.  The  latter  time  series  required  adoption  of  appropriate  hydrometeor  cloud 
models  which  are  described  in  Section  5. 

The  software  for  reading  the  SSM/I  tapes  was  implemented  on  the  Air  Force 
Geophysics  Laboratory  Interactive  Meteorological  System  (AIMS)  cluster  (Gustafson  and 
Felde,  1988).  The  unmodified  software  was  run  and  its  performance  verified.  Certain 
modifications  of  the  software  were  made  to  enhance  this  project's  performance.  The  first 
aspect  of  the  software  modification  concerned  the  performance  time  of  the  entire  package. 
In  its  unmodified  form,  the  tape  reading  package  took  about  8  hours  to  read  an  entire 
SSM/I  data  tape.  Since  there  are  32  tapes,  at  best  a  minimum  of  32  tape  reads  will  have  to 
be  performed,  which  would  translate  into  over  6  weeks  of  work  just  to  extract  the  data 
necessary  for  this  project.  After  closely  examining  the  individual  modules  and 
command/data  flow  design  of  the  software,  we  optimized  the  performance  time  such  that  it 
now  taken  slightly  more  than  1  hour  to  read  a  data  tape.  The  data  extracted  by  the  modified 
software  has  been  verified  and  validated  compared  to  data  extracted  by  the  unmodified 
package.  Software  was  developed  that:  (a)  loaded  the  site  locations  into  the  program,  and 
(b)  filtered  the  input  data  stream  from  any  SSM/I  data  tape  such  that  it  would  determine 
"local"  data  points  relevant  for  each  site,  and  group  them  into  respective  data  sets.  The  data 
is  binned  using  the  exact  latitude/longitude  coordinates  of  each  site  such  that  a  region 
defined  by  a  square,  whose  center  is  the  point  in  question,  is  generated.  Any  data  points 
which  fall  within  this  defined  region  are  considered  associated  with  the  point  in  question. 
Upon  completion  of  this  filtering,  a  file  for  each  region  is  written,  containing  the  data 
observed  by  the  SSM/I  for  the  time  period  in  question.  Software  to  read  the  SSM/I  tapes 
and  do  the  geographical  binning  by  site  is  contained  in  Appendix  A.  A  sample  output  from 
the  tape  reading  algorithm  is  contained  in  Appendix  B.  Output  are  the  time  (number  of 
seconds  since  the  SSM/I  sensor  was  first  turned  on),  latitude,  longitude  and  brightness 
temperatures  for  19.35  v,h,  22.235,  and  37  v,h  GHz.  The  85.5  GHz  channel  is  read  in  a 
second  pass  due  to  the  different  sampling. 

A  data  display  software  package  was  also  prepared  that  utilizes  NCAR  and  GKS 
routines  to  project  the  user's  choice  of  either  an  orthographic,  or  cylindrically  equidistant 
projection  of  the  earth,  and  then  superimposes  either  the  data  observed  for  a  specific  site,  or 
the  surface  scan  track  for  the  SSM/I  sensor  during  the  periods  in  question.  This  software 
is  given  in  Appendix  C.  Sample  map  output  illustrating  the  data  density  as  a  function  of 
site  at  various  scales  is  shown  in  Appendix  D. 

Software  was  also  written  to  evaluate  the  desired  climatologies.  These  elements 
include:  (a)  the  calculation  of  average  rain  rates  and  spatial  standard  deviations  within  each 
site  region,  and  (b)  the  evaluation  of  liquid  water  content  and  integrated  liquid  water  content 
time  series  based  on  the  adopted  hydrometeor  cloud  parameterization.  This  software  is  also 
provided  in  Appendix  A. 

3.2  SSM/I  Data  Set 

A  list  of  the  SSM/I  data  used  in  this  study  is  given  in  Appendix  E.  Provided  are: 

(a)  the  tape  number,  (b)  the  beginning  and  end  time  of  data  set,  and  (c)  the  number  of  files. 

The  data  files  were  acquired  as  part  of  this  study  and  are  available  for  further 
analysis.  The  software  capabilities  assembled  and  developed  during  the  course  of  this  pilot 
study,  i.e.  data  reading,  data  categorization,  data  plotting,  and  data  analysis,  should  be 
generally  applicable  to  the  analysis  of  this  SSM/I  data  set. 

While  no  comprehensive  examination  of  the  data  set  was  made  for  purposes  other 
than  those  of  this  study,  it  should  be  noted  that  this  global  data  set  should  be  extremely 
useful  for  a  variety  of  other  study  purposes. 
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4.  SSM/I  RAIN  RATE  CLIMATOLOGIES 

4. 1  Eurasian  Sites 

The  study  focused  on  eleven  USSR  stations  of  the  Environmental  Definition 
Program.  Data  on  these  sites  was  obtained  from  GL.  The  eleven  sites  are:  (a) 

Aktyubinsk,  (b)  Blagoveschensk,  (c)  Chita,  (d)  Kiev,  (e)  Leningrad,  (0  Mor  'ow,  (g) 
Murmansk,  (h)  Perm,  (i)  Semipalatinsk,  (j)  Simferopol,  and  (k)  Tashkent.  Tne  latitude 
and  longitude  coordinates  for  these  sites  used  in  the  study  can  be  found  in  subroutine 
"estreg"  in  Appendix  A. 

Of  these  sites,  Murmansk  is  the  northernmost  and  Tashkent  is  the  southernmost. 
Blagoveschensk  and  Chita  are  in  the  far  east  bordering  China.  Leningrad,  Moscow,  and 
Kiev  arc  the  most  western  sites.  Murmansk  is  the  only  coastal  site.  For  the  purposes  of 
this  study  all  of  the  sites  were  treated  as  land  based. 

This  study  did  not  call  for  the  collection  or  investigation  of  conventional  data 
sources  such  as  surface  and  upper  air  data  which  might  be  available  for  these  locations  or 
for  the  comparison  of  the  SSM/I  derived  rainfall  rate  climatologies  with  these  data.  In 
retrospect,  this  is  an  important  consideration.  Essentially,  the  subsequent  analysis  assumes 
that  the  statistical  relationship  between  rain  rate  and  brightness  temperature  inherent  in  the 
SSM/I  algorithm  regression  coefficients  (which  were  validated  over  the  United  States  and 
the  United  Kingdom),  are  equally  valid  over  the  set  of  Eurasian  sites. 

4.2  Site  Specific  Rain  Rate  Climatologies 

Referring  to  the  data  flow  diagram  presented  in  Figure  2,  the  next  step  in  the  data 
processing  is  the  coding  and  application  of  the  University  of  Wisconsin  SSM/I  rain  rate 
retrieval  algorithm  to  the  four  month  data  set.  This  has  been  accomplished.  Data  for  each 
of  the  eleven  sites  was  binned  according  to  lat/lo.ng  and  time  tag  and  rain  rate  retrievals 
were  performed  as  a  function  of  field-of-view  within  400  km  boxes  centered  at  each  site. 
This  bin  size  was  selected  somewhat  arbitrarily,  however,  consideration  was  given  to 
capturing  precipitation  events  of  synoptic  scale  which  passed  in  the  vicinity  of  the  site  as 
well  as  mesoscale/convective  activity.  Considerations  of  time-space  sampling  for  area- 
averaged  precipitation  (WMO,  1985)  were  considered  in  formulating  our  approach, 
however,  there  was  insufficient  time  to  fully  explore  these  issues. 

In  addition  to  the  spectral  information  available  in  the  SSM/I  brightness  temperature 
data  used  to  derive  the  surface  rainfall  rates,  it  was  recognized  that  the  spatial  distribution  of 
rain  within  the  binned  area  could  be  used  for  the  purpose  of  helping  to  characterize  the 
meteorological  properties  of  the  situation.  This  information  is  particularly  useful  to  aid  in 
the  selection  of  an  appropriate  parameterization  of  precipitation  liquid  water  content  (both 
vertical  distribution  and  integrated)  based  on  the  SSM/I  derived  surface  rain  rates.  It  is  the 
liquid  water  content  which  can  be  related  to  reentry  vehicle  erosion. 

For  this  reason,  both  the  average  daily  rain  rate  (defined  as  the  simple  arithmetic 
average  of  the  individual  SSM/I  derived  rain  rates  falling  within  the  site  specific  bin)  and 
the  spatial  standard  deviation  (SD)  within  the  bin  were  evaluated  to  analyze  the  spatial 
coherence  of  the  rain  rate  field.  These  data  were  plotted  as  time  series  to  produce  rain  rate 
climatologies  for  each  site  for  the  four  month  period.  Climatologies  for  each  site  are 
illustrated  in  Figures  3-13.  Illustrated  are  average  rain  rate  (Figs  3a-13a)  and  standard 
deviation  (Figs.  3b- 13b),  respectively.  The  time  series  are. labelled  in  days  from  1  June 
1989,  the  beginning  of  the  four  month  data  set.  All  time  series  plots  have  been  put  on  a 
common  scale  (truncated  at  2.5  mm/h)  so  that  intercomparisons  can  be  made. 

Examining  Figure  7a  for  Leningrad,  it  can  be  seen  that  there  are  obvious  high  (in  a 
relative  sense  averaged  over  400  km  squares)  rain  rate  situations  (e.g.  days  60-65)  and  low 
rain  rate  days  (e.g.  35-45).  The  standard  deviations  are  also  illustrative  (Fig.  7b).  For 
days  60-65,  the  high  rain  rates  are  accompanied  by  a  large  spatial  standard  deviation  (also 
days  27-30, 75,  and  82).  This  might  be  indicative  of  cellular  convection.  Moderate  rain 
rates  with  smaller  standard  deviations  might  indicate  more  uniform  precipitation.  The 
Leningrad  data  set  does  not  show  this  behavior. 
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Figure  4a.  Sample  time  series  afSSMI I  derived  rain  rate:  average  rain  rate. 


0  20  40  60  80  100  120 

Number  of  Ooys  from  01-JUN-1989 


Figure  4b.  Sample  time  series  ofSSM/I  derived  rain  rate:  standard  deviation. 
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Figure  7b.  Sample  time  series  qfSSM/I  derived  rain  rate:  standard  deviation. 
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Figure  9a.  Sample  time  series  qfSSMII  derived  rain  rate:  average  rain  rate. 
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Figure  9b.  Sample  time  series  ofSSM/1  derived  rain  rate:  standard  deviation. 
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Figure  10b.  Sample  time  series  of  SSM/I  derived  rain  rate:  standard  deviation. 
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Figure  lib.  Sample  time  series  ofSSM/I  derived  rain  rate:  standard  deviation. 
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Figure  12a.  Sample  time  series  cfSSMII  derived  rain  rate:  average  rain  rate. 


Figure  12b.  Sample  time  series  ofSSM/I  derived  rain  rate:  standard  deviation. 
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Figure  13b.  Sample  time  series  ofSSM/I  derived  rain  rate:  standard  deviation. 
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In  addition  to  the  time  series,  we  can  look  at  the  contoured  daily  SSM/I  derived  rain 
rates  for  given  sites.  This  gives  the  opportunity  to  study  the  structure  of  the  precipitation 
for  a  desired  day  at  the  resolution  of  the  SSM/I  footprints..  An  example  is  shown  in  Figure 
14  for  Moscow.  Based  on  the  time  series  data  given  in  Figures  8a, b,  the  figure  illustrates: 
(a)  a  low  rain,  high  SD  case  (day  15,  [41),  (b)  a  high  rain,  high  SD  case  (day  44, 10.41), 

(c)  a  high  rain,  low  SD  case  (day  61,  J0.2J),  and(d)  a  low  rain,  low  SD  case  (day  95, 

1.0.1 )).  The  numbers  in  brackets  are  the  contour  intervals  on  each  plot  in  mm/h.  A  heavy 
convective  cell  can  be  seen  in  Figure  14a,  but  there  is  little  precipitation  elsewhere  resulting 
in  a  veiy  low  average  rate  for  the  region.  Cells  of  moderate  intensity  can  be  seen  in  Figures 
14b, c  also.  In  these  cases  most  of  the  region  is  active  resulting  in  high  average  rainfall 
rates.  In  Figure  14d  by  comparison,  there  is  light  rain  throughout  the  region  (note  the 
contour  interval  changes).  TTie  importance  of  examining  both  the  average  rainfall  rate  and 
the  spatial  standard  deviation  in  characterizing  an  event  is  obvious. 

4.3  Discussion 

The  type  of  meteorological  event  resulting  in  precipitation  is  likely  to  be  related  to 
the  distribution  of  hydrometeor  signatures  within  the  region  examined.  High  standard 
deviation  indicates  that  the  precipitation  is  scattered  over  the  area  of  the  view  box,  e.g. 
convective  type  precipitation,  and  the  low  standard  deviation  means  the  rainfall  is  uniform 
over  the  area  of  the  view  box,  e.g.  stratiform  type  precipitation.  If  high  rainfall  rate  is 
associated  with  high  standard  deviation,  it  might  be  precipitation  from  individual 
thunderstorms.  The  passage  of  a  frontal  rainband,  on  the  other  hand,  could  be 
characterized  by  high  rainfall  rate  associated  with  low  standard  deviation.  The  association 
of  low  rainfall  rate  with  low  standard  deviation  could  mean  the  passage  of  a  warm  front, 
and/or  precipitation  from  stratiform  cloud.  If  low  rainfall  rate  associates  with  high  standard 
deviation,  it  could  mean  that  the  weather  system  which  causes  precipitation  passes  through 
only  part  of  the  view  box.  The  other  possibilities  of  low  rainfall  rate  with  high  standard 
deviation  could  indicate  that  the  weather  system  is  too  weak  to  produce  a  sufficient  amount 
of  precipitation  and/or  the  environment  is  too  dry,  thus  only  part  of  the  rainfall  could  be 
detected  at  the  surface  when  the  weather  system  passes  through  the  view  box. 

Regions  of  convective  and  stratiform  precipitation  are  hard  to  define,  because 
several  well-developed  and  decaying  convective  cells  might  still  overhang  into  the 
stratiform  region  (Tao  and  Simpson,  1989).  Therefore,  it  is  hard  to  use  the  rain  rate  to 
clearly  define  the  type  of  precipitation.  From  the  data,  we  can  see  that  the  two  stations 
located  at  the  eastern  part  of  the  USSR  have  higher  rainfall  rate  and  higher  standard 
deviation.  The  difference  between  these  two  stations  might  be  that  individual 
thunderstorms  occur  in  Blagoveschensk  more  than  Chita.  Although  both  data  show  that 
the  frontal  rainfall  occurs  very  often  at  both  sites.  It  might  be  that  the  Blagoveschensk  site 
is  much  closer  to  the  ocean.  Other  sites  have  less  rainfall  rate  and  show  the  situation  of 
steady  precipitation  occurs  more  frequently  in  these  locations,  despite  the  fact  that  some 
individual  thunderstorms  (maybe  afternoon  thunderstorms)  occurred  in  the  period. 

We  employ  these  spatial  coherence  concepts  to  the  application  of  our  precipitation 
cloud  models  (Section  5)  to  the  determination  of  integrated  hydrometeor  liquid  water 
content  in  Section  7. 

5 .  CONVECTIVE  CLOUD  MODEL  PARAMETERIZATION 

Task  2  required  definition  of  cloud  models  to  relate  SSM/I  derived  surface  rain  rates 
to  precipitation  liquid  water  content.  We  have  completed  this  task  by  defining 
parameterizations  of  liquid  water  content  vertical  distribution  and  total  integrated  liquid 
water  content  of  precipitation  for  stratiform  and  convective  rain  with  the  capability  to  decide 
between  the  two  using  the  spatial  coherence  data  derived  from  the  SSM/I  rain  rate  fields. 
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Figure  14.  Sample  contour  map  ofSSM/I  derived  rain  rate 
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5.1  Cloud  Process  Overv  iew 

The  formation  of  raindrops  or  ice  crystals  involves  complicated  microphysical  and 
dynamical  processes.  The  thermodynamical  effects  from  the  phase  changes  of  water  in 
turn  affect  the  evolution  of  the  weather  system  in  which  those  processes  are  embedded. 
Several  papers  have  studied  the  dependence  of  the  fonnation  of  precipitation  on 
microphysical  and  dynamical  processes.  The  results  show  that  both  the  microphysical  and 
the  dynamical  processes  are  important.  For  example,  the  growth  and  fallout  of  rain  can 
interact  with  the  updraft,  or  can  evaporate  in  the  subsaturated  downdraft  region,  which  in 
turn,  affects  formation  of  the  hydrometeors.  Since  measurement  of  the  microphysical 
processes  is  difficult,  the  only  way  to  simulate  these  processes  is  to  parameterize  them  with 
respect  to  the  large  scale  dynamical  motion.  Due  to  recent  theoretical  developments  and 
laboratory  experiments,  several  models  have  been  proposed  to  treat  these  processes  and 
some  of  the  results  have  been  compared  with  the  measurements  from  the  field  experiments. 

In  this  study  we  conducted  a  survey  of  different  hydrometeor  cloud  models 
focusing  on  those  which  can  be  used  to  simulate  the  vertical  distribution  of  liquid  water 
content  based  on  the  surface  precipitation  rate.  A  simple  parameterized  relationship  based 
on  climatological  statistics  is  derived  and  the  vertical  distribution  of  the  liquid  water  content 
is  expressed  as  a  function  of  height  and  surface  rainfall  rate. 


5 . 2  Review  of  Cloud  Vertical  Distribution  Models 
5.2. 1  Simulation  Studies 

A  high  correlation  between  cloud  top  height  and  rainfall  rate  has  been  shown  by 
Zawadzki  and  Ro  (1978)  and  Dennis  et  al.  (1975).  Adler  and  Mack  (1984)  used  a  one 
dimensional  cloud  model  to  study  the  relationship  of  the  thunderstonn  cloud  height-rainfall 
rate  and  the  cloud  height-volume  rainfall  rate  with  satellite  infrared  data.  Hie  variation  of 
the  vertical  velocity  and  precipitation  efficient  j  has  been  shown  to  dominate  both  the  slopes 
and  the  difference  of  the  two  curves  which  represent  the  relation  between  the  cloud  top 
height  and  the  rainfall  rate.  Information  on  the  convection  from  numerical  model  output  is 
needed  to  derive  the  cloud-rain  relationship,  and  the  vertical  shear  is  considered  to  be  the 
important  parameter  to  estimate  the  volume  rainfall  rate, 

Larger  scale  factors  such  as  the  synoptic  environment,  topographical  situation,  mill 
the  location  of  the  station  with  respect  to  the  thunderstorm  or  frontal  rainband  are  important 
in  detennining  the  desired  relationship.  Therefore,  empirical  rain  estimation  techniques 
developed  in  one  area  cannot  lie  applied  directly  to  other  areas.  Simple  adjustments  may  be 
inadequate  because  of  differences  in  slopes  of  rainfall  rate-height  in  different  locations.  To 
determine  convective  rain  rate  and  the  volume  rain  rate,  the  moisture  source,  the  vertical 
velocity,  and  the  rain  efficiency  are  important.  Additionally,  the  updraft  area  is  important  in 
determining  the  volume  rain  rate.  Ideally  all  of  these  data  should  be  used  to  adjust  the  value 
of  the  input  parameters  to  get  a  more  representative  profile.  Thus,  full  use  of  a  numerical 
cloud  model  with  the  inclusion  of  the  adequate  dynamical  processes  is  necessary. 

Kessler  (1959, 1961,  1963)devised  parameterized  equations  for  microphysical 
processes  (cloud  water  and  rain)  with  an  assumed  vertical  motion  profile  and  water 
generation  function.  Simpson  and  Wiggert  (1969)  used  a  one  dimensional  cloud  model 
with  the  Kessler  (1965)  type  of  parameterization  to  study  the  precipitation  in  tropical 
cumulus  clouds.  The  microphysical  processes  include  the  autoconversion  from  cloud 
water  to  precipitation  water,  collection  and  coalescence,  terminal  velocity,  and  fallout  of  the 
precipitation  with  the  evaporation  due  to  entrainment  of  drier  air  in  a  downdraft.  The  initial 
conditions  include  the  information  from  the  sounding  data:  the  saturation  at  cloud  base  or 
lifting  condensation  level  at  environment  temperature,  the  excess  of  the  temperature  at  cloud 
base,  and  the  vertical  velocity  at  the  cloud  base.  The  assumptions  made  with  the  above 
model  rule  out  the  feedback  between  the  microphysicai  and  dynamical  processes  in  the 
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model.  The  new  model  corrects  this  problem  and  the  results  are  more  reasonable. 

However,  the  treatment  of  the  ice  phase  is  insufficient  and,  therefore,  the  parameterization 
needs  to  be  updated  in  order  to  include  the  ice  phase  change  effect.  The  results  from  this 
study  reveal  that  more  complicated  microphysical  processes  should  be  introduced  and  the 
interaction  with  dynamical  effects  should  be  included.  The  amount  of  the  rain  which 
reaches  the  ground  as  precipitation  cannot  be  calculated  in  the  context  of  this  study  and, 
therefore,  we  cannot  use  the  information  of  the  surface  rainfall  rate  to  derive  the  vertical 
distribution  of  the  liquid  water  content. 

Cotton  (1972a,  1972b)  studied  precipitation  processes  within  the  supercooled 
cumuli  environment,  and  the  interaction  between  the  microphysical  process  tuid  cloud 
dynamics.  Precipitation  formation  in  warm  clouds  (1972a)  and  a  model  which  includes  the 
ice  phase  (1972b)  have  been  studied.  A  more  complicated  parameterization  to  simulate  the 
microphysical  processes  in  midlatitudes,  wliich  includes  the  processes  of  phase  change 
between  cloud  ice,  snow  and  graupel,  cloud  water  and  rain  water  has  been  developed  by 
Lin  et  al.  (1983).  This  parameterization  is  used  in  the  study  by  Rutledge  and  Hobbs  to 
investigate  the  precipitation  within  wann  clouds  (1983)  and  seeded  ice  phase  (1984) 
environments.  Three  different  precipitation  zones  categorized  according  to  the  surface 
precipitation  rate  «ire  described  in  the  study.  Their  results  are  in  a  good  agreement  with 
those  from  field  experiments. 

Tao  and  Simpson  (1989)  use  the  Kessler  type  of  microphysical  (cloud  water  and 
rain  water)  and  Lin  et  al.  (1983)  parameterizations  (cloud  water,  rain,  cloud  ice,  snow  and 
graupel)  to  study  the  structure  of  tropical  squall-type  convective  lines.  Two-dimensional 
models  and  three  dimensional  models  have  been  used.  The  role  of  the  ice-phase  and  the 
mesoscale  ascent  in  middle  and  high-levels  has  also  been  investigated.  The  model  output 
of  the  ice  runs  have  been  compared  with  that  of  the  ice-free  runs.  The  results  show  that  the 
ice-phase  microphysical  processes  are  crucial  for  a  realistic  stratiform  structure  and  its 
precipitation  statistics.  Also,  their  results  show  that  the  mesoscale  ascent  in  the  middle 
level  was  the  main  mechanism  responsible  for  the  extended  region  of  the  stratiform 
precipitation  at  the  rear  of  the  squall  line,  which  has  been  suggested  by  Rutledge  (1986). 

5.2.2  Measurement  Stud  ies 

■  Wei  et  al.  (1989)  use  five  different  methods  to  estimate  path-integrated  (columnar) 
cloud  liquid  water.  The  methods  include  one-channel  (31.65  GHz)  and  two-channel  (20.6 
GHz  and  31.65  GHz)  physical  retrievals,  the  standard  method  of  linear  statistical  inversion 
using  two  channels,  and  two  statistical  methods  that  proceed  from  an  initial  detennination 
of  several  empirical  regressions  between  measured  and  computed  quantities.  With 
brightness  temperature  data  and/or  the  absorption  coefficients  (for  oxygen,  the  water  vapor 
and  the  cloud),  the  optical  thickness  of  the  clear  air  and  the  cloud  water  can  be  calculated. 
The  calculation  of  cloud  water  content  with  microwave  radiometer  data  is  encouraging, 
however,  their  results  lack  comparison  with  independent  observations,  e.g.  radar. 

During  the  period  of  interest,  there  is  no  precipitation  occurring  in  the  study  of  Wei 
et  al.  (1989),  thus  they  only  estimate  the  cloud  water  content  without  estimating  the 
precipitation  amount.  The  study  needs  sounding  data  (vertical  distribution  of  pressure, 
temperature,  humidity,  and  cloud  water  content),  and  knowledge  of  the  absorption 
coefficients  to  resolve  the  answer.  Also,  a  hypothetical  liquid  water  profile  from  an 
archival  sounding  needs  to  be  specified  in  order  to  calculate  the  radiative  transfer.  These 
conditions  require  almost  the  same  effort  of  running  a  more  sophisticated  dynamical  cloud 
model. 

Heymsfield  and  Fulton  (1988)  studied  the  measurement  of  precipitation  in 
thunderstorms  from  high  altitude  remote  aircraft.  In  their  studies,  the  comparisons  between 
the  microwave  data  at  92  GHz  and  183  GHz  with  the  data  from  radar,  lidar  or  GOES  IR 
data  reveal  that  using  92  GHz  microwave  data  to  detect  the  rainfall  area  is  advantageous. 
The  liquid  water  content  (LWC)  or  ice  water  content  (IWC)  can  be  calculated  based  on  the 
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convective  area,  IWC  for  the  snow  aggregates  in  the  convective  region  and  IWC  in 
nonprecipitating  anvil  region.  The  relation  of  the  rainfall  rate  with  the  vertical  liquid  and/or 
ice  structure  is  unclear  in  this  study. 

Wilheit  et  al.  (1982)  used  the  19.35  GHz  and  183  GHz  data  from  a  microwave 
radiometer  to  study  the  precipitation  in  a  tropical  storm.  The  tendency  of  the  rainfall  rate 
can  be  matched  by  the  tendency  of  the  brightness  temperature  at  92  GHz,  and  the  183  GHz 
providing  some  information  on  the  vertical  extent  of  the  frozen  hydrometeors.  This  is  a 
more  direct  way  to  determine  the  vertical  distribution  of  the  liquid  water,  and  it  is  worth 
further  study. 

5.3  Discussion 

Simpson  and  Wiggert  (1969),  Cotton  (1972a, b),  Rutiedge  and  Hobbs  (1983, 

1984)  and  Tao  and  Simpson  (1989)  used  cloud  models  to  simulate  the  spatial  distribution 
of  the  liquid  water  distribution.  Since  the  vertical  distribution  of  liquid  water  depends  on 
the  stage  of  evolution,  the  position  away  from  the  core  of  the  updraft  and  the  type  of  the 
system,  the  dynamical  processes  are  as  important  as  the  microphysical  process.  Forvell 
and  Ogura  (1988)  found  that  the  addition  of  ice  was  responsible  for  the  achievement  of 
more  realistic  scale  features  in  the  convective  region  of  the  simulated  squall  line.  More 
complex  microphysical  process  which  included  the  ice-phase  parameterization  have  been 
simulated  by  Lin  et  al.  (1983).  Rutledge  and  Hobbs  (1984)  and  Tao  and  Simpson  (1989) 
show  that  the  model  output  is  in  good  agreement  with  the  field  experiment  data.  The  two- 
dimensional  cloud  model  of  Tao  and  Simpson  (1989)  is  capable  of  simulating  stratiform 
precipitation  and  the  areal  coverage  of  the  stratiform  region.  Therefore,  if  a  numerical 
approach  is  desired,  it  is  recommended  to  use  this  cloud  model  for  the  stratiform 
precipitation  cases. 

The  model  of  the  Rutledge  and  Hobbs  (1984)  is  good  in  the  situation  of  frontal  rain 
bands  and  convective  precipitation.  In  order  to  obtain  the  vertical  structure  of  the  liquid 
water  from  model  output,  several  initial  conditions  are  needed.  Using  this  model  requires 
the  specification  of  additional  sources  of  information.  Since  the  liquid  water  field  from 
model  output  is  specified  at  each  grid  point,  a  simple  relationship  Detween  the  vertical 
distribution  of  liquid  water  field  and  surface  rain  fall  rate  is  unnecessary.  If  we  have  more 
specific  description  of  the  liquid  water  field,  the  snow  field  and  the  cloud  water  field,  the 
total  liquid  water  content  in  a  column  can  be  obtained  by  integrating  vertically.  Also,  the 
rainfall  rate  is  determined  by  the  conditions  of  the  environment,  such  as  the  relative 
humidity,  the  wind  shear  and  updraft.  Therefore,  it  is  hard  to  find  a  simple  relationship 
which  will  satisfy  all  of  the  dynamical  conditions  in  determining  the  vertical  distribution  of 
liquid  water  depending  on  the  surface  rainfall  rate  alone. 

An  alternative  approach  to  derive  a  simple  relation  between  the  surface  rainfall  rate 
and  the  vertical  distribution  of  liquid  water  content  is  using  the  climatological  record.  From 
the  result  of  Adler  and  Mack  (1984),  the  height-rainfall  rate  relation  is  different  from  place 
to  place.  For  example,  high  rainfall  rate  for  moderate  and  small  storms  occurs  in  the 
coastal  regimes,  while  fairly  deep  convection  in  Midwest  thunderstorms  produces  only 
moderate  rain  rate  (Adler  and  Mack,  1984).  Therefore,  our  derivation  of  the  relationships 
between  the  surface  rainfall  rate  and  the  vertical  distribution  of  the  liquid  water  content 
based  on  the  climatological  data  depend  on  the  geographical  distribution  too. 

In  the  study  of  Adler  and  Mac  '1984),  the  statistical  retrievals  of  columnar  liquid 
water  and  water  vapor  are  found  to  be  more  accurate  than  physical  retrievals.  Therefore, 
we  will  use  the  climatological  data  to  derive  the  simple  relationship  based  on  the  vertical 
distribution  of  the  liquid  water  content  from  Falcone  et  al.  (1979).  The  important  features 
of  these  distributions  are  (see  Figure  15): 


(1)  A  maximum  liquid  water  content  is  located  between  the  cloud  top  and  the 
surface  (Rutledge  and  Hobbs,  1983, 1984;  Tao  and  Simpson,  1989).  From  the  model 
output,  the  height  of  the  maximum  liquid  water  content  is  strongly  influenced  by  the 
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vertical  velocity,  therefore,  the  maximum  height  of  the  cumulus  convection  is  normally 
higher  than  the  stratiform  precipitation,  and 

(2)  The  amount  of  the  liquid  water  content  is  nearly  uniform  below  the  cloud  base, 
or  the  liquid  water  content  below  the  freezing  level  in  the  convective  region  is  almost 
constant  with  height  (Tao  and  Simpson,  1989). 


5.4  Liquid  Water  Content/Rainrate  Relationships 

Based  on  the  results  from  the  diagnostic  studies,  a  statistical  polynomial  fit  like  the 
one  used  in  Wei  et  al.  (1989)  is  used  to  simulate  the  climatological  profiles  for  the  four 
precipitation  models  given  in  Falcone  et  al.  (1979).  Models  I-IV  (pp.  46-47)  correspond  to 
two  stratiform  (I,  II)  and  two  convective  (III,  IV)  models,  respectively.  The  following 
conditions  are  used  to  derive  the  coefficients  of  the  polynomial  fits: 


•  the  liquid  water  content  reaching  the  surface,  Ms,  is  derived  from  the  rain  rate 
retrieved  from  the  SSM/I  microwave, 


•  the  maximum  liquid  water  content  is  Mm  and  is  at  height  Zm, 


•  the  liquid  water  content  at  the  cloud  top  (Zt)  is  zero, 


•  the  first  derivation  of  the  polynomial  at  the  height  (Zm)  equals  zero,  and 


•  the  first  derivative  of  the  polynomial  at  the  surface  equals  zero. 


These  values  are  read  from  Figures  14-17  in  Falcone  et  al.  (1979),  see  Table  3. 

The  relationship  between  the  rainfall  rate  and  the  liquid  water  content 
can  be  found  in  Falcone  et  al.  (1979)  and  Simpson  and  Wiggert  (1969).  We  use  the 
relationship  of  Falcone  et  al.  (1979) 

M  =  0.07  *  RR  **0.83  , 

where  M  is  the  liquid  water  content  (g/m^)  and  RR  is  the  rainfall  rate  (mm/h),  and  apply  it 
at  the  surface  and  use  Ms  to  denote  the  liquid  water  content  here. 

A  fourth  order  polynomial  fit  is  used  to  simulate  convective  precipitation  (models 
III  and  IV).  The  result,  as  shown  in  Fig.  15a,  shows  the  simulation  of  the  vertical 
distribution  of  the  liquid  water  content  with  the  polynomial  equation  and  the  input  values  of 
Zt,  Zm,  Mm,  Ms  corresponding  to  the  four  precipitation  categories  from  Falcone  et  al. 
(1979)  (Models  I-IV,  pp  46-47),  given  in  Table  3.  The  cloud  top  of  curve  2  is  lower  than 
that  specified  by  Falcone  et  al.  (1979)  in  their  Fig.  15.  If  we  look  at  the  curves  carefully, 
we  can  find  that  the  slope  is  nearly  the  same  above  and  below  the  maximum  height.  In 
steady  rain  cases  like  stratiform  precipitation,  most  of  the  cloud  water  and  the  precipitation 
is  located  in  the  lower  troposphere.  The  maximum  liquid  water  accumulation  is  near  1  to  2 
km,  the  height  of  the  maximum  liquid  water  field  is  lower.  Thus,  if  the  slopes  above  and 
below  the  maximum  height  are  the  same,  the  liquid  water  content  will  vanish  at  a  height 
below  the  specified  cloud  top.  For  this  reason,  only  curves  3  and  4  are  used  to  model 
convective  rain. 

A  third  ord  :r  polynomial  equation  is  used  to  simulate  the  stratiform  precipitation. 
Three  tests  have  been  run  with  the  equations  satisfying  four  out  of  the  five  given 
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Table  3.  Polynomial  Fit  Data 


Curve 

RR 

(m/h) 

zt 

(km) 

Zm 

(km) 

(g/m3) 

1 

1.25 

1.5 

0.75 

1.1 

2 

5.00 

4.0 

1.50 

2.2 

3 

12.50 

5.0 

2.50 

4.8 

4 

15.00 

6.1 

3.50 

6.2 

conditions.  Test  one  satisfies  the  conditions  1,  3, 4, 5,  test  two  satisfies  the  condition  1, 2, 
3, 4  and  test  three  satisfies  conditions  1, 2, 3, 5.  The  test  two  results  (Figure  15b)  are  the 
best  fit  for  the  simulation  of  stratiform  precipitation. 

The  equation  for  the  vertical  distribution  for  the  convective  type  of  precipitation 
(e.g.  Falcone  models  III  and  IV),  is: 


LWC  =  a*z4  +  b*z3+c*z2  +  d. 


where 


,  1  f 

( 3znz ,2  -  2z,3  -  zj)ms  -  ( 3zmz 2  -  2z,3)mm] 

"  zmW-2zmW  +  zn3z,4l 

b  -  \ 

(4zm2z,2  - 2z,A  - 2zmA)ms  - (4zm3z,2  - 2z4)mn\ 

zmsz2-2zjz3  +  zm3z,^ 

-1 

[(zj  ’  4zmZ?  +  3z<)ms  +  {4zm  z3  -3z;-)mm\ 

zjz^-lzjz^zjz,4 

d  =  m. 

Curves  3  and  4  in  Figure  15a  provide  the  vertical  distribution  for  convective  situations. 
The  equation  for  the  stratiform  type  of  precipitation  (e.g.  Falcone  models  I  and  II) 


LWC  =  e*z3  +  f*z2  +  g*z  +  h, 


where 
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e  = 


=  Z.;2,(V  z ,f  '  Z-)i'".  +  Z'(2Z--2')m-] 


/  =  ,  ,'j - -  Z,3  ‘  22.‘)m,  -  (-^A  -  Z, ’)">„] 

ZmZt[Zr~Zm) 


8  =  - 7^ — -  3w,2  +  2z/)m,  +  (3zmz-  -  2z,3)mm] 

zm2,\z,‘zm) 


h  =  ms 

Curves  1  and  2  in  Figure  15b  provide  the  vertical  distribution  for  stratiform  situations  for 
different  input  parameters. 

The  vertical  integrated  liquid  water  content  is: 


ILWC  =  a/5  *  Zt5  +  b/4  *  Zt4  +  c/3  *  Zt3  +  e  *  Zt  +  Ms 


for  the  convective  precipitation,  and 

1LWC  =  e/4  *  Zt4  +  f/3  *  Zt3  +  g/2  *  Zt2  +  h  *  Zt  +  Ms 


for  the  stratiform  precipitation. 

Some  studies  show  that  the  height  of  the  cloud  top  and  the  surface  rainfall  rate  are 
closely  related  (Dennis  et  al.,  1975).  Thus  we  can  derive  the  statistics  of  the  relationship 
for  different  sites,  and  rewrite  the  cloud  top  as  a  function  of  the  surface  rainfall  rate  to 
derived  a  more  precise  vertical  distribution  of  the  liquid  water  field. 


5.5  Summary 

These  studies  reveal  that  the  formation  and  evolution  of  hydrometers  have  close 
relationships  with  both  microphysical  and  the  dynamical  processes.  The  inclusion  of  tire 
ice-phase  parameterization  is  important  in  simulating  these  processes. 

Since  the  vertical  distribution  of  the  liquid  and/or  ice  content  is  highly  dependent  on 
the  type  of  the  precipitation,  the  evolution  stage  of  the  thunderstorm,  and  the  distance  away 
from  the  core  of  the  thunderstorm  (in  the  area  near  the  convection  core  or  in  the  anvil  area), 
a  simple  relationship  between  the  surface  rainfall  rate  and  vertical  liquid  water  distribution 
is  quite  site  and  time  specific.  A  climatological  record  of  the  different  sites  should  be 
collected  and  used  to  derive  the  simulation  equation.  The  instant  observation  afforded  by 
the  satellite  qan  be  used  to  adjust  some  input  parameters,  e.g.  cloud  (op  and  surface  rainfall 
rate,  which  allow  the  simulation  to  lie  as  close  as  possible  to  the  real  situation.  t 

Hie  specific  simulated  liquid  water  field  cim  be  obtained  from  the  results  of  a 
numerical  cloud  model,  however,  this  requires  a  significant  amount  of  additional  data.  For 
this  reason,  we  have  adopted  vertical  profile  models  based  on  climatological  models 
requiring  a  minimum  of  input  data.  Tbe  relationships  provide  a  parameterization  of  verticiil 
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hydron.eteor  liquid  water  content  with  surface  rain  rate  and  other  input  sucli  as  the  height  of 
the  cloud  and  location  and  magnitude  of  the  maximum.  The  latter  are  also  provided  by  the 
Falcone  et  al.  (1979)  model  classes. 

6.  CLOUI)  LIQUID  WATER 

Recognizing  the  difficulties  associated  with  low  correlations  between  SSM/I 
brightness  temperatures  and  cloud  liquid  water  content  over  land  (and  snow),  the  SSM/I 
calibration  validation  team  recommended  that  retrievals  of  cloud  liquid  water  not  be 
attempted  over  land  (AlLshouse  et  al.,  1988).  In  this  study  we  have  specifically  undertaken 
to  explore  the  parameterization  of  hydrometeor  liquid  water  content  vertical  distribution  and 
integrated  amount  for  the  purpose  of  characterizing  intense  convective  activity.  In  the 
absence  of  other  data  sources  (e.g.  visible  and  infrared  satellite  data),  some  correlation  can 
be  inferred  between  precipitation  events  and  cloud  presence.  To  that  extent,  climatological 
cloud  models  can  be  applied  in  an  attempt  to  characterize  the  cloud  liquid  water  content 
based  on  SSM/I  detennination  of  the  hydrometeor  properties. 

To  this  end  we  have  adopted  the  cloud  models  proposed  by  Falcone  et  al.  (1979) 
which  accompany  his  respective  precipitation  event  vertical  profiles.  A  relationship 
between  rainfall  rate  and  cloud  hquid  water  of  the  fonn: 

M  =  0.05  *  RR  , 

is  proposed  to  provide  a  climatological  cloud  liquid  water  aqjunct  to  the  hydrometeor  liquid 
water  content  discussed  in  the  previous  section.  Following  Falcone  et  al  (1979),  the 
vertical  distribution  of  the  cloud  is  assumed  to  have  the  same  vertical  dependence  as  the 
rain,  except  that  its  magnitude  is  sca'sd  by  the  above  relationship.  Results  for  two  sites, 
Tashkent  and  Perm  are  illustrated  in  Figures  16  and  17. 

7 .  APPLICATION  TO  SELECTED  PRECIPITATION  EVENTS 

In  order  to  apply  the  liquid  water  content  parameterizations  for  convective  and 
stratiform  situations  defined  in  the  previous  sections,  criteria  had  to  be  developed  based  on 
the  mean  areal  rainfall  rates  and  spatial  standard  deviation  data  available  in  the  time  series. 
This  was  accomplished  by  examining  the  ensemble  properties  of  these  two  parameters  for 
the  four  month  study  period  using  a  simple  clustering  approach.  The  mean  areal  rainfall 
rates  and  spatial  standard  deviation  for  each  site  were  plotted  against  one  another.  Results 
are  given  in  Figures  18-28.  These  figures  were  plotted  on  the  same  scale  so  that 
comparisons  can  be  made  among  the  sites.  Note  that  area  averaging  reduces  local  values 
considerably.  Based  on  these  cluster  diagrams  an  approach  was  developed  to  differentiate 
between  stratiform  and  convective  precipitation  regimes. 

Notably  there  is  a  general  similarity  in  the  behavior  of  the  cluster  diagrams.  Two 
regions  are  definable  from  the  data.  In  the  first  region,  there  is  a  systemmatic  increase  of 
areal  standard  deviation  with  average  rain  rate.  For  some  sites  this  increase  is  gradual  turd 
near  linear  (see  Chita,  Figure  20),  while  for  other  sites  it  is  rapid  and  nonlinear  (see 
Leningrad,  Figure  22).  The  first  region  occupies  the  lower  left  quadrant  of  the  available 
data.  The  second  region  occurs  at  higher  rain  rates  and  standard  deviations  than  the  first 
region,  occupying  more  of  the  upper  right  quadrant  of  the  data.  The  break  points  defining 
the  transitions  from  region  one  to  region  two  vary  with  site.  Based  on  our  examination  of 
daily  contours,  we  have  defined  region  one  to  consist  largely  of  stratiform  events  and 
region  two  to  consist  largely  of  convective  events.  There  was  not  sufficient  time  in  the 
study  to  provide  detailed  analyses  to  support  this  hypothesis.  This  should  be  explored 
using  conventional  data  sources. 

The  implications  of  the  categorization  defined  above  based  on  examination  of  the 
cluster  diagrams,  is  to  provide  criteria  based  on  the  rain  rates  and  standard  deviations  to 
define  membership  in  one  of  the  clusters  (i.e.  either  stratiform  or  convective)  and  then  use 
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Figure  17.  Cloud  liquid  water  time  series  for  Perm. 
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Figure  19.  Cluster  diagram  for  Blagoveschensk 
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Figure  23.  Cluster  diagram  for  Moscow 
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Figure  26.  Cluster  diagram  for  Semipalatinsk 
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Figure  27.  Cluster  diagram  for  Simferopol 
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Figure  28.  Cluster  diagram  for  Tashkent 


the  appropriate  liquid  water  content  relationship  developed  in  Section  5  to  evaluate  time 
series  of  this  and  related  quantities  for  each  site. 

We  have  applied  this  approach  to  the  data  for  Moscow  displayed  in  Figure  23.  The 
Moscow  cluster  diagram  shows  a  region  of  low  standard  deviation  days  (SD  <  ~1.3)  with 
increasing  rain  rate  (region  1)  bounded  by  high  SD  days.  Thediighest  SD  days  are  those 
for  which  the  rain  rates  are  the  highest.  We  have  defined  this  latter  cluster  as  region  2  in 
the  context  of  the  previous  discussion.  Therefore,  the  criteria  adopted  for  this  site  is  that 
days  for  which  the  SD  is  less  than  1.3  are  defined  as  stratifonn  and  others  are  convective. 
The  time  series  corresponding  to  application  of  this  criterion  for  Moscow  is  shown  in 
Figure  29a.  For  comparison,  Figures  29b  and  29c  show  the  same  time  series  when  it  is 
assumed  that  all  precipitation  events  are  stratifonn  or  convective,  respectively. 

It  should  be  noted  that  this  criteria  applies  only  to  the  Moscow  time  series  and  that  a 
more  general  criteria  for  other  sites  has  not  been  developed.  Our  approach  would  be  to 
develop  the  generalized  criteria  based  on  nonnalized  values  using  the  mean  average  rainfall 
rate  (i.e.  the  mean  of  the  time  series  of  areal  averages)  and  mean  standard  deviation  (i.e.  the 
mean  of  the  time  series  of  spatial  standard  deviations). 


36 


Integrated  Liquid  Woter  Content  (g/m*^) 


13200 


13100 


13000 


12900 


12800  I _ I - 1 - 1 - 1 - — I L- 

0  20  40  60  80  100 

Number  of  Ooys  from  Ol-JUN-1989 


Figure  29c.  Time  series  ofILWC  for  Moscow  assuming  convective  only 


To  illustrate  application  of  this  approach  to  two  additional  sites  the  cluster  diagrams 
were  examined  and  Blagoveschensk  and  Leningrad  (Figures  19  and  22,  respectively)  were 
selected.  The  cluster  structure  for  Blagoveschensk  is  much  more  diffuse  than  for  Moscow. 
There  is  a  much  greater  variation  in  standard  deviation  which  is  qualitatively  manifested  in  a 
much  less  organized  stratifonn  arch.  In  this  case  we  have  based  the  definition  of 
convective  events  on  high  average  rainfall  rate  (>  1 .5  nun/li)  rather  than  on  standard 
deviation  as  was  done  in  the  Moscow  case.  The  results  of  applying  this  criterion  for  the 
determination  of  ILWC  for  the  Blagoveschensk  time  series  data  are  shown  in  Figures  30a- 
30c. 

The  cluster  data  for  Leningrad  exhibits  yet  another  form  of  behavior.  Rather  than 
being  diffuse  as  in  the  previous  case,  there  is  a  very  tight  cluster  of  low  rain  rate/  low 
standard  deviation  near  the  origin  with  some  outlying  points  of  higher  rain  rate  and 
standard  deviation.  We  have  defined  the  former  class  as  stratifonn  events  and  the  latter  as 
convective.  The  results  of  applying  these  criteria  for  the  determination  of  ILWC  for  the 
Leningrad  time  series  data  are  shown  in  Figures  31a-31c. 

8.  CONCLUSIONS 

This  study  has  established  the  feasibility  of  establishing  climatologies  of  stratifonn 
and  convective  precipitation  events  for  site  specific  regions  of  interests  based  on  the  use  of 
SSM/I  microwave  imager  brightness  temperature  data.  The  application  of  SSM/1 
microwave  imager  brightness  temperature  data  to  the  detennination  of  precipitation 
climatologies  for  eleven  selected  sites  has  been  investigated.  The  SSM/I  precipitation 
retrieval  algorithm  has  been  applied  to  the  detennination  of  surface  rainfall  rates  within  a 
400  km  region  surrounding  each  site  and  additionally  algorithms  have  been  developed  to 
provide  areal  averaged  rain  rates  and  spatial  standard  deviations.  A  novel  application  of  the 
spatial  standard  deviations  provides  infonnation  on  spatial  inhomogeneities  of  the  retrieved 
rain  rates  in  addition  to  exploiting  the  spectral  infonnation  content  of  the  microwave 
brightness  temperature  data. 
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Figure  30c.  Time  series  ofILWC  for  Blagoveschensk  assuming  convective  only 
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Figure  31b.  Time  series  of  ILWC  for  Leningrad  assuming  stratiform  only 


Figure  31c.  Time  series  of  ILWC  for  Leningrad  assuming  convective  only 


Time  series  of  the  areal  averaged  rain  rates  «uid  spati.il  standard  deviations  have 
been  calculated  for  a  four  month  summer  period  for  each  of  the  sites.  This  climatology  of 
precipitation  events  can  be  used  to  provide  an  indication  of  intense  convective  activity 
which  is  of  interest  in  the  determination  of  reentry  vehicle  erosion.  In  order  to  distinguish 
convective  precipitation  events,  the  areal  averaged  rainfall  rates  and  standard  deviations  for 
each  site  have  been  displayed  in  cluster  diagrams  for  the  entire  four  month  period.  Based 
on  these  cluster  diagrams,  a  set  of  preliminary  criteria  has  been  discussed  to  distinguish 
between  stratiform  and  convective  situatioas.  Parameterizations  of  the  vertical  distributions 
of  hydrometeor  liquid  water  content  for  both  stratiform  and  convective  situations  have  been 
developed  based  on  climatological  precipitation  cloud  models.  These  parameterizations 
were  developed  after  review  of  the  relevant  dynamical  cloud  model  literature.  Furthennore 
these  vertical  profiles  of  hydrometeor  liquid  water  content  have  been  integrated  to  provide 
functional  relationships  between  SSM/I  retrieved  surface  rain  rates  and  total  hydrometcor 
liquid  water  content.  Time  series  of  the  total  hydrometeor  liquid  water  content  have  also 
been  provided  for  selected  sites. 

It  is  recognized  that  cloud  liquid  water  content  profiles  are  also  desired  to  calculate 
reentry  vehicle  erosion  parameters.  While  cloud  liquid  water  content  is  not  directly 
retrievable  over  land  based  on  the  use  of  the  SSM/I  brightness  temperature  data,  model 
based  and  climatological  correlations  between  the  vertical  distribution  of  hydrometeor  liquid 
water  content  and  that  of  cloud  liquid  water  content  allows  for  a  nidimentary 
characterization  of  cloud  liquid  water.  A  simple  approach  correlating  hydrometeor  and 
cloud  liquid  water  vertical  distributions  has  been  explored  based  on  climatological  models 
available  from  the  microwave  attenuation  literature.  These  models  assume  that  the  cloud 
liquid  water  content  is  proportional  to  the  rainfall  rate  mid  that  the  shape  of  the  vertical 
distribution  of  cloud  liquid  water  is  similar  to  that  of  the  precipitation.  Using  this 
approach,  time  series  of  total  cloud  liquid  water  content  is  also  provided  for  selected  sites. 

9.  RECOMMENDATIONS 

As  proposed  this  study  has  focused  exclusively  on  the  use  of  SSM/I  microwave 
imager  brightness  temperature  data  to  develop  rain  rate  climatologies.  Precipitation  events 
are  indicative  of  meteorological  situations  which  are  inherently  problematic  with  respect  to 
reentry  vehicle  erosion.  A  novel  aspect  of  the  investigation  is  the  application  of  spatial 
coherence  concepts  to  the  characterization  of  precipitation  regimes  such  as  stratiform  and 
convective  which  determine  the  vertical  distribution  of  hydrometeor  liquid  water  content. 
Rain  also  acts  as  a  surrogate  for  die  presence  of  cloud  which  is  not  easily  measured  over 
land  with  microwave  sensors. 

A  number  of  avenues  of  additional  study  are  recommended  based  on  the  results 
reported  here: 

•  More  work  should  be  done  to  exploit  the  identification  of  convective  precipitation 
using  the  spatial  infonnation  content  of  the  microwave  data.  The  examination  of  the  spatial 
coherence  properties  of  each  site  for  die  study  period  should  be  augmented  with 
conventional  observations  to  help  in  stratifying  the  observed  arches  and  clusters  according 
to  the  nature  of  die  precipitation  event.  For  example,  surface  and  upper  air  data  should  be 
available  from  ETAC  to  be  used  in  the  analysis; 

•  As  noted  in  Section  5  it  is  possible  to  identify  numerical  models  of  precipitation 
dynamics  to  suppoil  parameterization  of  rain  rate/liquid  water  content  relationships.  In  this 
study  we  resorted  to  climatological  relationships  for  expediency.  A  fully  interactive 
mesoscale  model  could  be  used  as  a  vehicle  for  satellite  data  fusion  and  retrieval  purposes; 

•  The  algorithms  developed  here  can  lie  readily  applied  to  the  examination  of  alternate 
sites  and  changes  in  the  statistical  computations  such  as  the  time  period  mid  spatial 
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averaging  field.  Additionally,  although  bandpass  limited  by  the  inherent  instantaneous 
fields-of-view  of  the  SSM/I  data,  the  spatial  power  spectrum  of  precipitation  events  could 
be  determined.  This  would  provide  additional  degrees  of  freedom  beyond  the  simple 
measure  provided  by  the  spatial  standard  deviation; 

•  The  clustering  of  rain  rate  and  spatial  standard  deviations  can  be  extended  to  use 
surface  observations  to  further  stratify  the  precipitation  regime  characterization.  For 
example,  frontal  and  various  thunderstonn  categories  could  be  added  along  with  their 
respective  vertical  distribution  models.  Available  surface  observations  should  also  be 
employed  to  verify  tuid  validate  the  SSM/1  time  series.  This  was  not  investigated  here; 

•  Data  fusion  with  sensor  data  sources  colocated  aboard  the  DMSP  spacecraft  (OLS, 
SSM/T,  SSM/T-2)  should  be  investigated.  Visible  and  infrared  imagery  provide  a  means 
to  specify  the  presence  of  cloud,  cloud  cover,  cloud  type,  and  cloud  top  height.  The 
evolution  of  multispectral  imagery  will  provide  the  capability  to  identify  high,  middle,  and 
low  cloud  in  conjunction  with  nephanaiysis  techniques  (d'Entremont,  1986).  The  use  of 
microwave  and  millimeter  wave  sounder  channels  provides  a  means  to  sense  the  vertical 
distribution  of  liquid  water  content  for  precipitation  analogous  to  the  "slicing"  approach  for 
clouds  employed  using  infrared  sounder  channels.  The  SSM/T  microwave  temperature 
sounder  provides  a  signature  in  each  channel  which  is  proportional  to  the  presence  of 
precipitation  elements  at  altitudes  characteristic  of  its  weighting  function.  Grasiewski  and 
Staelin  (1989)  have  suggested  utilizing  the  impact  of  precipitation  on  microwave  oxygen 
band  (60GHz)  and  line  (1 18GHz)  absorption  to  obtain  precipitation  vertical  structure 
information.  An  application  of  these  concepts  to  the  profiling  of  hydrometeor  liquid  water 
from  DMSP  SSM/T  data  could  be  fruitful  in  this  regard.  This  would  give  a  direct 
measurement  which  frees  one  from  the  climatological  precipitation  vertical  distribution 
models  employed  in  this  study.  The  vertical  distribution  of  cloud  liquid  water  can  then  be 
obtained  using  the  simple  climatological  based  correlations  suggested  in  this  study; 

•  The  SSM/T-2  moisture  sounder  is  also  sensitive  to  precipitation,  however,  due  to 
its  shorter  wavelength  response,  it  also  provides  a  cloud  liquid  water  signature  (Isaacs  and 
Deblonde,  1987).  Although  the  number  of  available  channels  (five)  is  not  as  extensive  as 
that  potentially  available  from  an  infrared  sounder,  the  millimeter  wave  channels  respond  to 
cloud  liquid  water  content  which  can  be  directly  related  to  cloud  vertical  structure  and 
erosion  parameter  indices.  The  problem  of  high,  variable  surface  emissivity  inherent  in  the 
use  of  the  SSM/I  imager  data  is  mitigated  by  the  use  of  the  SSM/T-2  sounder  data.  This  is 
due  to  the  vertical  resolution  properties  of  the  sounding  weighting  functions  which 
selectively  sense  the  middle  and  upper  tropospheric  cloud  liquid  water  (for  increasingly 
absorbing  channels  of  the  sounder); 

•  Finally,  it  is  recommended  that  the  use  of  data  available  from  civilian  satellites 
including  both  the  Tiros  polar  orbiters  and  the  GOES  geosynchronous  platforms  be 
investigated.  Use  of  multispectral  cloud  imagery  from  the  Tiros  Advanced  Very  High 
Resolution  Radiometer  (AVHRR)  can  be  used  for  multispectral  nephanaiysis.  High 
resolution  Infrared  Sounder  (HIRS)  data  can  be  used  for  applying  the  C02  slicing  method. 
The  GOES  VAS  (VISSR  Atmospheric  sounder)  provides  visible  and  infrared  imagery  and 
infrared  sounder  data.  GOES  infrared  data  could  be  particularly  useful.  The  specific 
approach  adopted  could  be  modeled  on  that  described  by  Adler  and  Negri,  1988).  They 
used  GOES  imager  data  to  delineate  convective  rain  areas  by  searching  for  minima  in  the 
equivalent  black  body  brightness  temperature  (EBBT)  field  and  then  assigned  rain  rates 
based  on  the  results  from  a  one  dimensional  cloud  model  which  provided  the  relationship 
between  convective  development  (cloud  top  height)  and  the  resulting  rain  rate.  Stratiform 
rain  rates  were  delineated  based  on  EBBT  tlireshold  criteria. 


43 


10.  ACKNOWLEDGEMENT 

We  thank  Dr.  Ken  Champion  and  G.  Felde  of  the  Geophysics  Laboratory  for  their 

assistance  during  this  effort. 

11.  REFERENCES 

Adler,  R.  F.  and  R.  A.  Mack,  1984:  Thunderstonn  cloud  height-rainfall  rate  relations  for 
use  with  satellite  rainfall  estimation  techniques.  J.  Climate  Appl.  Meteor..  23,  280-296. 

Adler,  R.  F.  and  A.  J.  Negri,  1988:  A  Satellite  Infrared  Technique  to  Estimate  Tropical 
Convective  and  Stratiform  Rainfall.,  A.  Appl.  Met..  22,  30-5 1 . 

Alishouse,  J.C.,  1983:  Total  precipitable  water  and  rainfall  detennination  from  the 
SEAS  AT  scanning  multichannel  microwave  radiometer,  J.  Geophvs.  Res..  88(C3), 
1929-1935. 

Alishouse,  J.C,  S.  Snyder,  J.  Vongsathom,  and  R.  Ferraro  (1988):  Water  vapor  and 
cloud  water  validation  in  DMSP  special  sensor  microwave/imager  calihration/validation. 
J.  Hollinger  (ed.).  Naval  Research  Laboratory  (Draft). 

Atlas,  D.  and  O.  Thiele  (eds)  ,1981:  Precipitation  measurements  from  space.  Workshop 
report,  Goddard  Space  Flight  Center,  Greenbelt,  MD. 

Bunting,  J.T.,  and  C.N.  Touart,  1980:  Horizontal  scale  variations  in  satellite  estimates  of 
weather  erosion  parameters  for  reentry  systems.  AFGL-TR-80-0209.  ADA092447 

Conover,  J.1L,  and  J.T.  Bunting,  1977:  Estimates  from  satellites  of  weather  erosion 
parameters  for  reentry  systems.  AFGL-TR -77-0260.  ADA053654 

Cotton,  W.R.,  1972a:  Numerical  simulation  of  precipitation  development  in  supercooled 
cumuli,  1.  Monthly  Weather  Rev..  100. 757-764. 

Cotton,  W.R.,  1972b:  Numerical  simulation  of  precipitation  development  in  supercooled 
cumuli,  2.  Monthly  Weather  Rev..  100. 764-784. 

Dennis,  A.  S.,  A.  Koscielski,  D.  E.  Cain,  J.  H.  Hirsch  and  P.  L.  Smith,  1975:  Analysis 
of  radar  observations  of  a  randomized  cloud  seeding  experiment.  J.  Appl.  Meteor..  14. 
897-908. 

d'Entremont,  R.P.,  1986:  Low-  and  midlevel  cloud  analysis  using  nighttime  multispectral 
imagery.  J.  Clim.  Appl.  Met..  25.  1853-1869. 

Falcone,  V.J.,  L.  ,V.  Abreu,  and  E.P.  Shettle,  1979:  Atmospheric  attenuation  of  millimeter 
and  submillimeter  waves:  Models  and  computer  code,  AFGL-TR-79-0253.  |NT1S  AD 
A084485J 

Feteris,  P.J.,  A.S.  Lisa,  and  A.J.  Bussey,  1976:  Environmental  definition  program  cross 
sectional  analysis:  Summary  of  data  and  analysis  techniques.  AFCRL-TR -76-0002. 
ADA021908 

Forvell,  R.  G.,  and  Y.  Ogura,  1988:  Numerical  simulation  of  a  midlatitude  squall  line  in 
two  dimensions.  J.  Atmos.  Sci..  45,  3846-3879. 


44 


Gloerson,  P.,  D.J.  Cavalieri,  A.T.C.  Chang,  T.T.  Wilheit,  W.J.  Campbell,  O.M. 
Johannessen,  K.B.  Katsaros,  K.F.  Kunzi,  D.B.  Ross,  D.  Staelin,  E.P.L.  Windsor, 
F.Y.  Barath,  P.  Gudmandsen,  c.  Langham,  and  R.O.Ramseier,  1984:  A  summary  of 
results  from  the  first  NIMBUS  7  SMMR  observations,  J.  Geophvs.  Res..  89(D4), 
5335-5344. 


Grasiewski,  A.J.  and  D.H.  Staelin  (1989):  Numerical  modeling  of  passive  microwave  02 
observations  over  precipitation.  Submitted  to  Radio  Science. 

Gustafson,  G.B.  and  G.W.  Felde,  1988:  Interactive  satellite  image  processing  applied  to 
cloud  detection.  Preprints:  4th  conference  on  Interactive  Information  Processing 
Systems.  AMS,  Boston,  MA. 

Hardy,  K.R.,  1979:  Meteorological  data  analysis  in  support  of  the  ABRES  program. 
AFGL-TR-79-0084.  ADA070946 

Heymsfield,  G.M.  «uid  R.  Fulton,  1988:  Comparison  of  high-altitude  remote  aircraft 
measurements  with  the  radar  stnicture  of  an  Oklahoma  thunder-storm:  Implications  for 
precipitation  estimation  from  space.  Mo.  Wea.  Rev..  1 16.  1 157. 

Huang,  R.,  and  K.N.  Liou,  1983:  Polarized  microwave  radiation  transfer  in  precipitating 
cloudy  atmospheres:  Applications  to  window  frequencies,  J.  Geophvs.  Res..  88(C6). 
3885-3893. 

Isaacs,  R.  G.  and  G.  Deblonde,  1987:  Millimeter  wave  moisture  sounding:  the  effect  of 
clouds.  Radio  Science.  22.  3,  367-377. 

Isaacs,  R.G.,  1989:  A  unified  retrieval  methodology  for  the  DMSP  meteorological 
sensors,  in  RSRM  '87:  Advances  in  Remote  Sensing  Retrieval  Methods.  A.  Dcepak, 
H.E.  Fleming,  and  J.S.  Theon  (eds.).  A.  Deepak  Publishing,  pp.  203-214. 

Isaacs,  R.G.,  Y.-Q.  Jin,  R.D.  Worsham,  G.  Deblonde,  and  V.J.  Falcone,  Jr.  (1989):  The 
RADTRAJN  microwave  surface  emission  models,  IEEE  Transactions  on  Geoscience  and 
Remote  Sensing.  27,  4,  pp.433-440. 

Isaacs,  R.G.,  R.N.  Hoffman,  and  L.D.  Kaplan  (1986):  Satellite  remote  sensing  of 
meteorological  parameters  for  global  numerical  weather  prediction,  Reviews  of 
Geophysics.  24.  4,  pp.  701-743. 

Jin,  Y.Q.  and  Isaacs,  R.G.,  1985:  Thermal  radiative  transfer  for  inhomogeneous 

nonisothermal  atmospheric  precipitation  and  statistical  retrieval  of  related  parameters.  In 
Remote  Sensing  Instrumentation:  Technology  for  Science  and  Applications. 

Proceedings  of  the  1985  International  Geoscience  and  Remote  Sensing  Symposium 
(IGARSS’85).  Amherst,  MA,  7-9  October,  1985.  Univ.  Massachusetts,  IEEE. 

Jin,  Y.Q.  and  R.G.  Isaacs,  1987:  Simulation  and  Statistical  Retrieval  for  Inhomogeneous 
Non-isothermal  Atmospheric  Precipitation.  J.  Quantitative  Spectroscopy  and  Radiative 
Transfer.  37.  5,  461-469. 

Kessler,  E.,  1959:  Kinematic  relations  between  wind  and  precipitation  distributions,  I.  T 
Meteor..  16,  630-637. 


45 


Kessler,  E.,  1961:  Kinematic  relations  between  wind  and  precipitation  distributions,  II.  L. 
Meteor.,  jl,  510-525. 

Kessler,  E.,  1963:  Elementary  theoiy  of  associations  between  atmospheric  motions  and 
distributions  of  water  content.  Mon.  Wea.  Rev..  91, 13-27. 

Kessler,  E.,  1965:  Microphysical  parameters  in  relation  to  tropical  cloud  and  precipitation 
distributions  and  their  modification.  GcofLsica  International.  79-88. 

Kessler,  E.,  1969:  On  the  distribution  and  continuity  of  water  substance  in  atmospheric 
circulation.  Meteor.  Monogr..  32,  Amer.  Meteor.  Soc.,  84  pp. 

Lin,  Y.  -L.,  R.  D.  Farley  and  H.  D.  Orville,  1983:  Bulk  parameterization  of  the  snow  field 
in  a  cloud  model.  J.  Climate  Appl.  Meteor..  22, 1065-1092. 

Liou,  K.N.,  G.C.  Aufderhaar,  and  P.T.  Nipko,  1981 :  Some  examples  of  the  effects  of 
clouds  and  precipitation  on  the  temperature  profile  retrieval  for  DMSP  SSM/T 
microwave  sounders,  J.  Appl.  Meteorol..  20,  821-825. 

Lipes,  R.G.,  1982:  Description  of  SEASAT  radiometer  status  and  results.  J.  Geophvs. 
Res..  87(C5),  3385-3396. 

Lo,  R.  C.,  1983:  A  Comprehensive  Description  of  the  Mission  Sensor  Microwave  Imager 
(SSM/1)  Environmental  Parameter  Extraction  Algorithm,  NRL  Memorandum  Report 
5199.  Naval  Research  Laboratory,  52  pp. 

Marshall,  J.S.  and  W.M.  Palmer,  1948:  'Hie  distribution  of  raindrops  with  size.  J. 
Meteor..  5.  165166. 

NASA,  1976:  Satellite  derive  global  oceanic  rainfall  atlas.  NASA  SP-410. 

Negri,  A.J.  and  R.F.  Adler,  1987a:  Infrared  and  Visible  Satellite  Rain  Estimation.  Part  I: 
A  Grid  Cell  Approach.  J.  Clim.  and  Appl.  Met..  26.1565-1576. 

Negri,  A.J.  and  R.F.  Adler,  1987b:  Infrared  and  Visible  Satellite  Rain  Estimation.  Part  II: 
A  Cloud  Definition  Approach.  J.  Clim.  and  Appl.  Opt..  26.  1565-1576. 

Olson,  W.S.,  F.J.  LaFontaine,  W.L.  Smith,  R.T.  Merrill,  B.A.  Roth,  and  'I'.  Achtor, 
(1988):  Recommended  algorithms  for  the  retrievjil  of  rainfall  rates  at  midlatitudes  using 
the  SSM/I  (DMSP-8).  In  DMSP  special  sensor  microwave/imaner 
calibration/validation.  J.  Hollinger  (ed.)  Naval  Research  Laboratory  (Draft). 

Peirce,  R.M.,  R.W.  Lenhard  and  B.D.  Weiss,  1975:  Comparison  study  of  models  used  to 
prescribe  hydrometeor  water  content  values.  AFCRL-TR-75-0470.  ADAO 19633 

Rutledge,  S.  A.  tuid  P.  V.  Hobbs,  1983:  The  mesoscale  and  microscale  structure  and 
organization  of  clouds  and  precipitation  in  midlatitude  cyclones.  VIII:  A  model  for  the 
"seeder-feeder"  process  in  warm-frontal  rainbands.  J.  Atmos.  ScL  40.  1 1 85- 1 206. 

Rutledge,  S.A.,  and  P.V.  Hobbs,  1984:  The  mesoscale  and  microscale  structure  and 
organization  of  clouds  and  precipitation  in  mid-latitude  cyclones.  XII:  A  diagnostic 
modeling  study  of  precipitation  development  in  narrower  frontal  rainbands.  J.  Atmos. 
Sci..  41.  2949-2972. 


46 


Rutledge,  S.  A.,  1986:  A  diagnostic  modeling  study  of  the  stratiform  region  associated 
with  a  tropical  squall  line.  J.  Atmos.  Sci..  43.  1356-1377. 

Savage,  R.  C.,  1981:  The  DMSP  Microwave  Imager  (SSM/I):  Unfinished  Plans.  In 
Precipitation  Measurements  from  Space.  D.  Atlas  and  O.  Theile,  (eds).  pp  246-252. 

Savage,  R.C.  and  J.A.  and  Weinman,  1975:  Preliminary  calculations  of  the  upwelling 
radiance  from  rain  clouds  at  37.0  and  19.35  GHz.  Bulletin  of  the  American 
Meteorological  Society.  56.  1272-1274. 

Savage,  R.  C.,  L.  A.  Burgess,  and  B.  W.  Shipley,  1987:  The  Defense  Meteorological 
Satellite  Program  (DMSP)  Microwave  Imager  (SSM/I)  in  Passive  Microwave  Observing 
from  Environmental  Satellites.  J.  C.  Fisher  (ed.).  NOAA  Technical  Report  NESDIS 

35.  pp.  186-202. 

Simpson,  J.,  and  V.  Wiggeit,  1969:  Models  of  precipitating  cumulus  towers.  Monthly 
Weather  Rev..  97.  471-489. 

Simpson,  J.,  R.  H.  Simpson,  D.  A.  Andrews  and  M.  A.  Eaton,  1965:  Experimental 
cumulus  dynamics.  Reviews  of  Geophysics,  3,  387-431. 

Spencer,  R.W.  and  D.A.  Santek,  1985:  Measuring  the  global  distribution  of  intense 
convection  over  land  with  passive  microwave  radiometry,  J.  Cli.  Met.. .24.  860. 

Tao,  W.-K.,  and  J.  Simpson,  1989:  Modeling  study  of  a  tropical  squall-type  convective 
line.  J.  Atmos.  Sci..  46.  177-202. 

Tsang,  L.,  and  J.A.  Kong,  1977:  Theory  for  thermal  microwave  emission  from  a  bounded 
medium  containing  spherical  scatterers.  J.  App!.  Physics.  4£,  3593-3599. 

Ulbrich,  C.W.,  1983:  Natural  variations  in  the  analytical  form  of  the  raindrop  size 
distribution.  J.  Appl.  Meteorol..  22.  1764-1775. 

Weaver,  R.,  C.  Morris,  R.  G.  Barry,  1987:  Passive  Microwave  Data  for  Snow  and  Ice 
Research:  Planned  Products  from  the  DMSP  SSM/I  System,  EOS.  29  September,  p. 
769. 

Wei.  Leighton,  mid  Rogers,  1989:  A  comparison  of  several  radiometric  methods  of 
deducing  path  integrated  cloud  liquid  water.  JAOT.  6,  1001-1012. 

Weinman,  J.A.,  and  P.J.  Guetter,  1977:  Detennination  of  rainfall  distribution  from 
microwave  radiation  measured  by  the  NIMBUS  6  ESMR,  J.  Appl.  Meteorol..  16. 437- 
442. 

Weinman,  J.A.,  and  T.T.  Wilheit,  1981.  A  survey  of  passive  microwave  and  hybrid 
remote  sensing  of  precipitation,  in  Precipitation  Measurements  from  Space  Workshop 
Report,  edited  by  D.  Atlas  mid  O.W.  Thiele,  GSFC,  pp.  D-173-D-175. 

Wilheit,  T.T.,  A.T.C.  Chang.  M.S.  V.  Rao.  E.B.  Rodgers,  and  J.S.  Theon,  1977: 
Satellite  technique  for  quantitatively  mapping  rainfall  rates  over  the  oceans,  J.  Appl. 
Meteorol..  16,  551-560. 


Wilheit,  T.T.,  A.T.C.  Chang,  M.S.V.  Rao,  E.B.  Rodgers,  R.A.  Nieman,  B.M.  Krupp, 
A.S.  Milman,  J.S.  Stratigos  mid  H.  Siddalingaiali,  1982:  Microwave  radiometric 


47 


observations  near  19.35, 92,  and  183  GHz  of  Piecipitation  in  Tropical  Storm  Cora.  J. 
Appl.  Met..  21.  1137. 

Willis,  P.T.,  and  P.  Tattleman,  1989:  Drop-size  distributions  associated  with  intense 
rainfall.  J.  Appl.  Meteor..  28.  3-15. 

WMO,  1985:  Review  of  requirements  for  area  averaged  precipitation  data,  etc.  WCP-100. 
WMO/TD-NO.  1 15,  World  Meteorological  Organization. 

Zawadzki,  1. 1.,  and  C.  V.  Ro,  1978:  Correlations  between  maximum  rate  of  precipitation 
and  mesoscale  parameters.  J.  Appl.  Meteor..  17,  1327-1344. 


48 


Appendix  A  SSM/I  Data  Extraction  Software 
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0001  SUBROUTINE  GTLAT  {  GTL ) 

0002 

0003  c..*#........*...*.......******^****.#.*******.****.****.*. 

0004  C  CALCULATES  LATITUDE  OF  GROUND  TRACK  POINT  FROM  REGRESSION 

0005  C  ON  SCAN  PIXELS  3,63  AND  125  (IE.  2,32,63  IN  64  PIX  SCAN) 
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Appendix  B  Sample  Output  from  the  Tape  Reading  Algorithm 
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Time 

Longitude 

Latitude 

Brightness  Temperature  (K) 

(Seconds) 

(Degrees) 

(Degrees) 

19.35v 

19.35h 

22.235 

37v 

37h 

86543433 

55.68 

38.46 

270.60 

266.80 

270.40 

269.50 

265.90 

86543433 

55.81 

38.79 

272.40 

269.30 

271.40 

269.70 

266.70 

86543433 

55.93 

39.13 

272.10 

269.30 

271.20 

270.50 

267.60 

86543437 

55.02 

37.48 

271.10 

268.80 

271.60 

269.40 

266.50 

86543437 

55.17 

37.78 

270.80 

268.20 

270.80 

269.30 

266.10 

86543437 

55.31 

38.08 

270.00 

265.10 

271.40 

268.40 

264.30 

86543437 

55.46 

38.39 

271.10 

264.50 

269.30 

267.40 

263.30 

86543437 

55.59 

38.72 

271.60 

267.50 

270.10 

269.60 

266.60 

86543437 

55.71 

39.06 

271.50 

268.60 

272.00 

269.30 

266.00 

86543440 

54.80 

37.42 

271.90 

268.70 

270.60 

268.60 

265.40 

86543440 

54.95 

37.72 

271.90 

268.60 

271.60 

269.30 

266.20 

86543440 

55.09 

38.02 

271.50 

267.60 

269.70 

268.30 

266.00 

86543440 

55.24 

38.33 

271.00 

264.90 

269.10 

267.60 

264.00 

86543440 

55.37 

38.66 

270.50 

265.60 

270.20 

267.90 

263.40 

86543440 

55.49 

39.00 

•  271.90 

268:60 

270.30 

269.10 

266.00 

86543444 

54.58 

37.36 

271.60 

268.20 

269.10 

268.10 

264.90 

86543444 

54.73 

37.66 

271.20 

268.20 

269.90 

268.50 

264.70 

86543444 

54.87 

37.96 

271.20 

267.70 

269.90 

268.50 

266.10 

86543444 

55.02 

38.26 

271.20 

265.90 

268.10 

268.10 

264.90 

86543444 

55.15 

38.59 

270.60 

265.20 

269.60 

267.80 

262.90 

86543444 

55.27 

38.93 

270.90 

266.80 

270..Q0 

268.80 

264.90 

86543448 

54.35 

37.30 

270.30 

266.70 

268.70 

267.60 

264.80 

86543448 

54.50 

37.60 

271.10 

266.90 

267.60 

267.90 

266.00 

86543448 

54.65 

37.90 

271.40 

267.80 

270.40 

268.80 

264.80 

86543448 

54.80 

38.20 

270.30 

267.30 

269.90 

267.60 

264.40 

86543448 

54.93 

38.53 

270.00 

264.90 

268.10 

267.10 

263.30 

86543448 

55.05 

38.86 

270.40 

264.70 

269.80 

266.90 

262.80 

86543448 

55.18 

39.19 

270.90 

266.00 

269.30 

268.30 

265.10 

86543452 

54.28 

37.54 

269.80 

266.70 

268.50 

267.20 

263.70 

86543452 

54,43 

37.83 

270.40 

266.20 

268.30 

268.10 

264.90 

86543452 

54.58 

38.13 

270.40 

267.20 

270.40 

268.40 

264.90 

86543452 

54.71 

38.46 

269.70 

264.90 

268. -30 

267.20 

262.70 

86543452 

54.83 

38.79 

268.80 

263.60 

268.70 

266.90 

262.00 

86543452 

54.96 

39.12 

269.70 

263.40 

268.50 

266.90 

261.70 

86543456 

54.21 

37.77 

270.10 

266.60 

268.40 

266.40 

262.80 

86543456 

54.36 

38.07 

270.30 

266.30 

269.90 

267.30 

264.40 

86543456 

54.49 

38.39 

270.20 

265.70 

270.10 

266.90 

264.40 

86543456 

54.61 

38.72 

269.10 

264 . 30 

268.20 

265.60 

262.10 

86543456 

54.74 

39.05 

269.20 

263.80 

269.00 

266.60 

261.80 

86543459 

54.27 

38.33 

270.20 

265.70 

268.80 

267.10 

263.00 

86543459 

54.39 

38.66 

269.30 

264.00 

268.40 

266.80 

263.50 

86543459 

54.52 

38.98 

269.00 

263.70 

268.00 

266.40 

262.50 

86543463 

54.17 

38.59 

269.60 

264.50 

267.40 

266.50 

262.00 

86543463 

54.30 

38.91 

269.70 

264.30 

267.60 

266.30 

262.50 

86543467 

54.20 

39.17 

269.80 

266.00 

268.10 

267.30 

264.20 

86720726 

57.73 

39.10 

267.70 

264.20 

264.90 

264.20 

262.80 

86720730 

57.70 

38.66 

266.20 

263.70 

267.30 

265.60 

264.30 

86720730 

57.52 

38.93 

267.50 

265.20 

267.50 

265.40 

263.20 

86720730 

57.34 

39.18 

267.20 

264.50 

266.80 

263.70 

261.30 

86720733 

57.68 

38.22 

264.70 

259.80 

264.00 

263.80 

260.50 

86720733 

57.50 

38.48 

267.60 

264.50 

267.30 

266.40 

265.20 

86720733 

57.32 

38.75 

268.30 

265.50 

267.30 

265.50 

264.00 

86720733 

57.14 

39.00 

267.80 

264.90 

266.00 

265.10 

263.20 

86720737 

57.63 

37.74 

263.50 

,  258.70 

264.50 

261.80 

259.70 

86720737 

57.47 

38.04 

263.50 

258.00 

264.90 

262.50 

257.20 

86720737 

57.29 

38.30 

267.50 

264.60 

265.70 

266.80 

264.90 

86720737 

57.11 

38.57 

268.30 

265.20 

267.40 

26b. 00 

265.10 

86720737 

56.93 

38.82 

267.10 

263.70 

266.50 

264.30 

261.20 
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Appendix  D  Sample  Mapping  Software  Output 
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Appendix  E  SSM/I  Data  Catalog 
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E-1 


SSM/I  Data  Received 


Tape  8 

Beginning  of  dataset 

End  of  dataset 

#  files 

01 

01-Jun-89 : 00 : 07 : 31 

05-Jun-89 : 00 : 57 : 56 

56 

02 

05-Jun-89 : 00 : 57 : 56 

09-Jun-89 : 00 : 06 : 24 

20 

03 

09-Jun-89 : 00 : 06 : 24 

13-Jun-89 : 00 : 56 : 52 

55 

04 

13-Jun-89 : 00 : 56 : 53 

16-Jun-89 : 00 : 18 : 21 

42 

05 

16-Jun-89 : 00 : 18 : 21 

20-Jun-89 : 01 : 08 : 51 

49 

06 

20-Jun-89 : 01 : 08 : 51 

24-Jun-89 : 00 : 17 : 17 

54 

07 

24-Jun-89 : 00 : 17  : 18 

28- Jun-89 : 01 : 07 : 33 

49 

08 

28-Jun-89 : 01 : 07 : 33 

01-Jul -89:00:28:44 

42 

09 

01-Jul -89:00:28:55 

05-Ju 1-89:01:18:59 

49 

10 

05-Jul-89 : 01 : 18 : 58 

09-Jul-89:00:27:20 

53 

11 

09-Jul-85:00:27:20 

13-Jul-89:01:17:35 

54 

12 

13-Jul-89 : 01 : 17 : 35 

16-Jul-89:00:38:48 

42 

13 

16-Jul-89 : 00 : 38 : 40 

20- Jul-89 : 01 : 29 : 01 

57 

14 

20-Jul-89:01:29:01 

24-Jul-89:00:37:17 

12 

15 

24-Jul-89:00:37:17 

28- Jul-89 : 01 : 27 : 89 

32 

16 

28-Jul-89:01:27:29 

Ol-Aug-89 : 00 : 35 : 45 

56 

17 

Ol-Aug-89 : 00  :  35  :  45 

05-Aug-89 : 01 : 25 :  54 

57 

18 

05-Aug-89 : 01 : 25 : 54 

09 -Aug-8 9 : 00 : 34 : 07 

51 

19 

09-Aug-89 : 00 : 34 : 07 

13 -Aug- 8 9 : 01 : 24 : 18 

57 

20 

13-Aug-89 : 01 : 24 : 19 

16 -Aug-8 9 : 00 : 45 : 28 

42 

21 

16-Aug-89 : 00 : 45 : 28 

20-Aug-89 : 01 : 35 : 36 

56 

22 

20-Aug-89 : 01 : 35 : 36 

2 4 -Aug-8 9 : 00 :  43  :  47 

53 

23 

24-Aug-89 : 00 : 4  3  :  47 

28-Aug-89 : 01 : 33 : 51 

57 

24 

28-Aug-89 : 01 : 33 : 51 

01-Sep- 89:00:41:48 

56 

25 

01-Sep- 89:00:41:46 

05-Sep-89:01:31:45 

55 

26 

05-Sep-89:01:31:45 

09-Sep- 89:00:39:47 

56 

77 

09-Sep-89:00:39:47 

13-Sep-89:01:29:43 

57 

28 

13-Sep-89:01:29:43 

16-Sep-89:00:50:41 

42 

29 

16-Sep-89 : 00 : 50 : 41 

20-Sep-89:01:40:33 

57 

30 

20-Sep-89:01:40:33 

24-Sep-89:00:48:26 

56 

31 

2 4 -Sep-8 9 : 00 : 48 : 26 

28-Sep-89:01:38:15 

57 

32 

28-Sep-89 : 01 : 38  : 15 

Ol-Oct-89 : 00 : 59 : 09 

40 

There  is  no  data  missing  from  our  requested  set. 
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